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A EXTRACT 

This  report  analyzes  and  reviews  the  processes  involved 
in  the  generation  of  modulations  in  short  wave  structure  by 
internal  waves,  ine  emphasis  is  on  field  conditions  where, 
in  general,  there  is  wind  blowing,  a mean  current  with  verti- 
cal structure,  long  and  short  waves  at  different  directions 
to  the  internal  wave  and  possibly  surface  films.  The  appli- 
cability cr  otherwise  of  theoretical  analyses  and  laboratory 
experiments  to  field  situations  is  discussed,  the  problem 
being  the  variety  of  interacting  phenomena  fodnd  under  natural 
conditions.  At  the  end  of  each  main  section  of  the  report  a 
summary  is  given  of  the  principal  conclusions  of  -chat  section. 


In  Section  2,  the  mean  distribution  in  the  vertical  of 
the  horizontal  current  is  discussed.  In  the  top  few  milli- 
meters there  is  a surface  wind  drift,  generally  about  3 per- 
cent of  the  wind  speed.  Below  this,  the  current  profile  is 
highly  variable  in  time  and  in  location,  but  is  essentially 
horizontal,  with  a much  greater  horizontal  length  scale  tho.n 
vertical  scale.  Response  times  are  estimated.  Section  3 
considers  the  variation  in  surface  current  produced  by  internal 
waves  in  the  presence  cf  such  a mean  current  distribution.  It 
is  shown  that,  for  a given  internal  wave,  the  variations  are 
amplified  by  a surface  current  in  the  same  direction  as  the 
internal  wave-  propagation,  reduced  if  opposed. 

The  formation  of  surface  films  in  the  presence  of  internal 
waves,  currents,  winds,  surface  -waves,  etc.,  is  discussed  in 
Section  h and  limits  are  given  for  conditions  under  which 
satura ted  films 


should  bo  found 
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Section  5 reviews  known  results  concerning  the  direct  in- 
teraction between  a surface  current  distribution  and  surface 
waves  and  considers  their  applicabilit5'  to  field  situations. 
Except  in  very  unusual  circumstances  (very  fast  internal  waves , 
surface  current  in  the  opposite  direction)  the  wave  lengths  in^ 
volved  in  direct  resonance  are  so  short  that  energy  input  from 
the  wind  will  m--.sk  this  process  — resonance  is  not  considered 
to  be  a significant  process  in  generating  short  wave  modulations, 
Longer,  faster  surface  waves,  whose  time  constant  for  energy - 
input  from  the  wind  is  much  greater,  also  experience  modulations' 
in  their  propagation  through  an  internal  wave  field,  though  to 
a considerably  smaller  extent.  However,  these  modulations  do 
influence  strongly  the  structure  that  short  waves  can  maintain 
against  breaking.  The  process  of  wave  breaking  is  considered 
in  Section  6 and  estimates  are  made  of  the  modulations  in  short 
wave  energy  density  resulting  from  this  indirect  interaction 
in  the  presence  of  wind.  The  results  are  summarized  in  Sec- 
tion 6 . . 
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1.  INTRODUCTION  AND  OUTLINE 

An  tinderstanding  of  the  processes  involved  in  the  produc- 
tion of  surface  effects  by  internal  gravity  waves  or  other  sub- 
surface disturbances  is  essential  for  reliable  prediction,  but 
this  understanding  lias  remained  elusive.  A variety  of  physical 
phenomena  are  potentially  involved  — non-resonant  coupling 
between  short  surface  waves  and  the  moving  disturbance  pattern; 
resonant  interactions  which  are  more  selective  in  their  occur- 
rence but,  under  idealized  conditions  at  lc-ast,  potentially 
stronger;  indirect  coupling  between  bhort -waves  and  the  moving 
pattern,  either  via  longer  'waves  or  swell  or  via  sub-surface 
wind-induced  currents;  modulations  in  the  pattern  or  intensity 
of  small  scale  wave  breaking  again  either  directly  or  enhanced 
by  an  intermediary  such  as  surface  wind  drift;  modulations  in 
the  density  of  surface  material  with  consequent  variations  in 
energy  loss  rate  and  energy  density  of  short  waves.  Under 
natural  conditions  some  or  all  of  these  may  occur  in  varying 
combinations;  what  is  demonstrable  in  the  laboratory  under 
carefully  controlled  conditions  and  justified  theoretically  by 
a correctly  argued  analysis  may  be  overwhelmed  in  the  field  by 
a quite  different  process  where  there  is  a wide  range  of  per- 
turbations of  many  kinds.  Certainly  some  of  the  early  observa- 
tions of  LaFond  (19o2)  on  the  occurrence  and  movements  of 
slicks  pointed  rather  definitely  to  an  association  with  in- 
ternal waves,  but  under  other  sets  of  natural  conditions,  other 
balances  obtain  and  the  associations  are  no  longer  clear. 
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The  aim  of  this  report  is  to  summarize  what  we  know  about 
the  processes  that  occur  near  the  sea  surface  that  may  be 
relevant  to  these  questions  and  to  evaluate  their  interplay 
and  their  combined  consequences  insofar  as  they  produce  an 
observable  variation  In  surface  properties  as  a result  of  a 
sub-surface  disturbance.  What  theoretical  models  do  we  have, 
that  offer  guidance  to  observation  in  the  laboratory  and  in 
the  field?  What  physical  processes  are  responsible  for  the 
phenomena  already  observed  in  the  laboratory?  Which  Of  these 
are  relevant  in  the  field  and  what  are  the  optimum  conditions 
for  detectability  in  the  field  — which  property  of  what  pro- 
vides modulations  under  the  widest  range  of  conditions  and 
what  is  the  best  way  of  measuring  it?  This  is  the  range  of 
questions  to  which  this  report  is  addressed  in  the  hope  that 
It  will  place  much  of  the  excellent  work  already  done  into 
context  and  indicate  directions  where  future  work  may  be  most 
fruitful.  ■' 

The  wind  is  generally  blowing  over  the  ocean.  The  sur- 
face drag  produces  a thin  surface  drift  layer  and,  if  it  has 
been  blowing  for  a substantial  time,  sub-surface  wind-induced 
currents  are  present  also.  Waves  exist  over  a continuous  range 
of  scales  that  can  extend  from  capillaries  to  long  swells;  in 
an  active  wind-generated  sea,  waves  are  breaking  sporadically, 
sometimes  entraining  air  and  forming  a visible  white  cap  but 
sometimes  on  such  a small  scale  that  no  air  entrainment  takes 
place.  Generally  evaporation  produces  a cooled  surface  ’skin1; 
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when  wave  breaking  occurs  this  skin  is  locally  and  abruptly 
disrupted  so  that  the  local  surface  temperature,  measured 
rad  iome  trie  all  -r,  is  also  disrupted  abruptly.  Surface  films, 
usually  of  organic  origin,  are  present  almost  ubiquitiously . 
Superimpose  on  all  this  a very  small  amplitude,  large-scale 
disturbance  generated  by  an  internal  wave  — what  are  the  con- 
sequent large  scale  modulations? 

These  phenomena  will  be  considered  in  turn  with  what  we 
know  about  their  self-  and  mutual  interactions.  First  is  the 
wind-induced  drift,  both  in  the  thin  surface  layer  and  on  a 
larger  scale  that  may  extend -to  depths  of  many  meters.  What 
are  the  response  times  of  these  to  changes  in  surface  stress? 
Secondly,  how  do  they  respond  to  variations  in  the  underlying 
current,  whether  steady  currents  or  the  moving  orbital ' current 
patterns  of  internal  waves?  Third  are  the  surface  films,  their 
properties  and  the  variations  in  film  density  induced  by  sur- 
face motion,  whether  produced  by  surface  waves  or  internal 
waves,  possibly  augmented  by  wind-drift  currents.  What  is  the 
back  reaction  on  the  waves  themselves,  either  to  variations  in 
the  damping  of  short  components  or  to  an  influence  of  the 
density  on  the  vigor  of  wave  breaking?  Fourth,  we  consider 
the  possible  modulations  in  surface  wave  energy  density  pro- 
duced by  near-surface  current  distributions,  either  by  resonance 
(with  maybe  ’blockage')  or  non-resonance.  What  influence  does 
the  wind  have  on  these  modulations?  What  effect  does  the  wind 
drift  in  the  thin  surface  layer  have  on  wave  propagation  and 
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the  occurrence  of  these  modulations?  Fifth,  what  about  wave 
breaking?  We  show  that  this  is  influenced  strongly  by  the  thin 
wind-induced  surface  layer  and  that  the  short  wave  energy  den- 
sity that  can  be  maintained  against  breaking  depends  on  the 
combined  influence  of  longer  waves  and  the  thin  surface  layer. 

In  all  of  these  questions,  the  aim  is  to  draw  together  the 
results  from  theory,  from  the  laboratory  and  from  the  field  to 
try  to  develop  a consistent  picture  or  to  pinpoint  where  in-  * . 
consistencies  lie.  The  emphasis  will  be  on  the  surface  effects 
themselves,  rather  than  on  the  mode  of  generation  of  the  sub- 
surface disturbance  that  produces  them.  We  will  assume  that, 
by  one  means  or  another,  we  have  a disturbance  at  a given  depth; 
what  are  its  surface  manifestations?  But  these  are  enough 
questions,  and  it  is  time  to  look  for  answers. 
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2.  MEAN  VELOCITY  STRUCTURE  ABOVE  THE  THERMOCLINE 
2.1  The  Surface  Drift 

When  the  v;ind  blows  over  the  surface  of  the  sea,  the  stress 
in  the  air  or  the  water  is  communicated  by  the  Reynolds  stresses 
associated  with  the  almost  invariably  turbulent  motion.  In  the 
absence  of  wave  breaking,  however,  the  tangential  stress  across 
the  interface  must  be  communicated  from  the  air  to  the  water  by 
molecular  viscosity,  since  tne  Reynolds  stress  at  the  interface 
vanishes.  This  requires  a substantial  velocity  gradient  at  the 
surface  in  both  the  air  and  the  water: 


where  xQ  is  the  tangential  stress  at  the  interface  and  p and  v 
are  the  density  and  kinematic  viscosity  of  either  the  air  or 
water.  With  increasing  distance  on  either  side  of  the  inter- 
face, the  Reynolds  stress  -puw  increases  rapidly  and  soon 
dominates  so  that  the  mean  velocity  gradient  reduces  rather 
abruptly.  This  viscous  sublayer  is  well  known  in  aerodynamics; 

i 

over  a smooth  plate  its  thickness  is  approximately  10v(p/to)2 
and  the  velocity  at  its  outer  edge  relative  to  the  plate  is 
about  10(r  /p)H  = 10u#,  where  u#  is  the  friction  velocity  of  the 
air  flew.  Even  a smooth  water  surface  is  mobile  in  the  sense 
that  tangential  motions  are  possible  so  that  the  build,  up  of 
the  Reynolds  stress  above  and  below  the  air -water  interface 
would  be  expected  to  be  more  rapid  and  the  viscous  sublayer 
consequently  thinner  than  over  a flat  plate.  Nevertheless  on 
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dimensional  grounds,  the  sealing  is  the  .samej  the  thickness 
of  this  layer  in  the  air  is  of  order  v&/u * and  the  velocity 
scale  is  u#,  where  v&  is  the  viscosity  of  the  air.  The  stress 
is  continuous  across  the  Interface  so  that  the  velocity  scale 
In  the  water  motion 


w 


# 


3.4  x 10‘2  u# 


[2.21 


where  pa  and  p represent  the  air  and  water  densities  respec- 
tively. The  thickness  of  the  viscous  sublayer  is  proportional 
to  v /w  . Relative  to  the  speed  of  the  interface  itself,  the 
velocity  profiles  in  the  air  and  the  water  have  an  anti-symmetry 
but  with  different  length  and  velocity  scales.  The  ratio  of 
the  length  scales  may  be  uncertain  if  the  flow  is  aerodynamically 
rough,  but  the  ratio  of  the  velocity  scales  in  water  to  air  is 
simply  wyut  = 3J1  percent.  Consequently,  the  speed  q of  the 
water  at  the  interface  relative  to  the  water  below  (where  the 
mean  velocity  profile  is  quite  flat)  is  predicted  to  be  about 

3 percent  of  the  air  speed  as  usually  measured  at,  say,  10 

meters  where  again  the  mean  profile  is  relatively  flat. 

In  summary,  then,  we  would  expect  on  these  grounds  a thin 
wind-drift  layer  occupying  the  top  f ew  millimeters  of  the  water 
surface  across  which  the  velocity  difference  is  of  the  order  3 
percent  of  the  wind  speed.  The  existence  of  this  layer  has  only 
a minor  influence  on  the  propagation  speed  of  waves  whose  wave- 
length Is  large  compared  with  the  layer  depth  (Fhillips,  19?3)j 
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since  the  dispersion  relation  is  the  result  of  an  integral 
balance  over  the  whole  wave  motion.  However,  as  will  be  seen 
later,  it  does  influence  very  strongly  the  condition  for  wave 
breaking  that  intimately  involves  the  motion  right  at  the  water 
surface.  If  the  waves  do  break,  the  viscous  sublayer  is  dis- 
rupted locally  and  must  be  re-established  by  the  surface  stress. 
The  time  to  accomplish  this  can  be  estimated  by  equating  the 
impulse  (force  times  time)  supplied  per  unit  area  by  the  surface 
wind  stress  to  the  momentum  (water  density  times  depth  times 
velocity)  that  has  to  be  acquired  by  the  layer.  Thus, 


w 

a"*  ~e  'w  w*  W* 


P u *T  * P* 


so  that  the  re-establishment  time 


P v 

w 


p U* 
ra  * 


[2.3] 


For  wind  speeds  of  10  knots  or  greater  this  characteristic  time 
is  very  short,  a fraction  of  a second,  and  only  in  very  light 
winds  (of  order  2 knots)  is  the  ro-establishment  time  of  the 
order  1 second  or  greater. 

Hov/ever,  for  small  changes  in  surface  stress  or  in  u#,  the 
response  time  for  readjustment  may  be  substantially  greater. 

If  u^  increases  to  ut  + bu^,  the  readjustment  time  is 
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or  greater  than  [2.3]  by  & factor  (uy£6u^).  Consequently  in 
a non-breaking  wave,  when  the  tangential  wind  stress  at  the 
crests  is  greater  by  a few  percent  (Brooke  Benjamin,  1959)  than 
in  the  troughs,  the  response  of  the  thin  surface  layer  to  the 
stress  variation  is  negligible  for  waves  of  one  second  period 
or  less,  provided  the  wind  speed  is  not  more  than  about  15  or 
20  knots.  On  the  other  hand.  Equation  [2.3]  shows  that  if  the 
sublayer  is  disrupted  entirely,  it  is  essentially  re-established 
very  quickly  by  the  wind  stress. 

Measurements  of  the  surface  wind  drift  in  the  field  were 
made  some  years  ago  by  Keulcgr.n  (1951) i ho  reported  a surface 
drift  of  approximately  3*3  percent  of  i he  wind  speed,  Wu  (1968) 
in  a careful  laboratory  study  measured  not  only  the  surface 
drift  but,  by  using  spherical  floats  of  different  sizes,  was 
able  to  confirm  the  sharp  sub-surface  velocity  gradient  and  its 
variation  with  wind  speed.  Figure  1 shows  his  measurements  of 
dU/dz  as  a function  of  Uq2,  the  square  of  the  centerline  wind 
speed  in  his  tunnel,  which  is  very  closely  proportional  to  the 
surface  stress.  Over  most  of  the  range  the  two  are  propor- 
tional., confirming  that  the  viscous  balance  of  Equation  [2.1] 
is  indeed  relevant.*  When,  at  the  highest  wind  speeds,  the 

Incidentally,  the  mean  viscous  stress  in  the  wind  drift  layer 

can  be  estimated  from  the  best  fit  line  of  Figure  1 to  be 

0.8  x 10“3  p U 2 in  the  absence  of  wave  breaking,  somewhat 
ra  o 

greater  than  Wu  estimates  for  the  total  stress  from  his  air 
flow  measurements.  Since  there  is  also  momentum  flux  to  the 
waves,  it  should  be  less  -*■'  the  reason  for  this  discrepancy  is 
unknown . 
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waves  were  breaking,  the  mean  velocity  gradient  was  reduced 
primarily  as  a result  of  the  sporadic  patches  of  vigorous 
turbulence  and  the  locally  increased  Reynolds  stresses.  Wu's 
results  for  the  surface  drift  are  shown  in  Figure  2,  together 
with  others  obtained  by  Phillips  and  Banner  (1974)  in  a dif- 
ferent laboratory.  These  observations  and  experiments  certainly 
confirm  the  existence  of  this  wind-drift  laydr,  both  in  the 
laboratory  ana  the  field,  the  magnitude  of  tip1 surface  drift  q 
being  about  3 to  4 percent  of  the  wind  speed 
theory  predicts. 


as  indeed  the  simple 


2.2  Sub-Surface  Currents 

The  wind  also  generates  sub -surface  currents  on  a larger 
scale  of  course.  The  vertical  gradient  of  Reynolds  stress  is 
balanced  by  the  horizontal  acceleration  of  the  fluid  and/or  by 
the  Coriolis  acceleration.  In  the  open  sea,  With  a gradually 
developing  'wind  ever  a large  area,  the  mean  horizontal  current 
distribution  can  approach  the  classical  Ekman  spiral,  with  a 
depth  scale  of  the  order  w^/fi  (perhaps  100  meters  or  so),  where 
n is  the  vertical  component  of  the  earth's  rotation.  The  con- 
ditions necessary  for  this,  however,  are  extreme  and  are  not 
generally  to  be  expected.  The  vertical  profile  of  mean  hori- 
zontal current  in  the  mixed  layer  over  depths  from  1 meter,  say, 
to  the  the rmoc line  certainly  depends  on  the  past  history  in 
space  and  time  of  winds  in  the  vicinity,  on  whether  or  not  the 
point  of  interest  is  near  a shoreline,  a current  system  or  a 
detached  eddy,  and  must  be  expected  to  be  highly  variable. 
Curiously  enough,  few  detailed  field  measurements  of  this  kind 
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soem  to  have  been  taken  except  for  a limited  number  in  par- 
ticular current  systems  like  the  Pacific  Equatorial  System 
(Montgomery  and  Stroup,  19v>P).  Experiments  such  as  MODE 
concentrate  on  greater  depths  and  have  sparse  coverage  in 
this  range.  Nevertheless,  it  is  probably  a safe  assertion 
that  there  is  no  reason  to  expect  that  the  vertical  profile  of 
horizontal  current  in  the  mixed  layer  is  either  zero,  uniform 
or  a classical  Ekman  spiral ; we  must  allow  for  considerable 
variability. 

Suppose,  then,  that  there  is  a current  distribution  in 
the  mixed  layer  in  which  the  Reynolds  stress  gradients  arc 
balanced  by  the  mean  acceleration  and  by  Coriolis  accelerations. 
Suppose  new  that  the  flow  field  is  disturbed  by  an  internal 
wave  or  other  perturbation  in  the  neighborhood  of  the  thermo- 
cline.  Both  the  distribution  of  Reynolds  stresses  and  of  mean 
current  will  be  perturbed;  if  the  perturbation  is  instantaneous, 
the  altered  Reynolds  stresses  have  no  time  to  react  back  on 
the  mean  flow,  but  the  perturbation  is  not  instantaneous  and 
the  question  arises:  Under  what  conditions  will  the  perturba- 

tion in  Reynolds  stress  influence  the  induced  perturbation  in 
the  current? 

The  question  can  be  posed  analytically  as  follows:  We 

take  a frame  of  reference  moving  with  the  internal  wave  train 
as  in  Figure  3 and  define  phase  averaged  quaiU  '.ties  < > as 
averages  along  the  y-direc tion , parallel  to  the  internal  wave 
crests.  The  mean  velocity  field  defined  in  this  way  is  U(x,z), 
V(x,z),  W(x,z)  and  the  total  velocity,  including  turbulent  and 
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surface  wave  fluctuations  is  U(x,z)  + u(x,y,z,t),  etc.  The 
upper  mixed  layer  above  the  thermocline  at  a depth  z = -d  is 
generally  cf  almost  uniform  density  and  the  time  scale  of  the 
internal  waves  of  concern  are  short  compared  with  one  pendulum 
day,  sc  that  in  this  region,  Coriolis  and  buoyancy  effects  can 
reasonably  be  neglected.  The  momentum  equations  in  the  x-  and 
z-directions  are: 


..  bU  bfJ  * > 

ox  be  ax  ^ ^ as  -u,;> 


v ■ <UW>  + <v;2> 

<?*  cx  az 


i a<o>  ) 
P bx  j 

i a<p> 

p bz 


[2*3] 


whence,  by  cross -differentiation 
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[2.7] 


is  the  mean  vorticity  in  the  flew.  The  right  hand  side  of 
[?.6]  represents  the  influence  on  the  flow  of  the  perturbations 
in  Reynolds  stresses  that  result  from  the  internal  wave  induced 
disturbance.  In  the  absence  of  the  internal  wave  but  with  the 
horizontal  current  structure,  tj  <±  U(z),  V = 0,  <uw>  = const  and 
<u^>,  <w2>  = functions  of  z ad  one;  all  terms  vanish  identically. 


/ 


-14- 


HYDRONAUTICS,  Incorporated 


The  Reynolds  stress  perturbations  cannot  be  estimated 
precisely  from  analytical  considerations.  It  has  been  argued 
in  the  context  of  air  flow  ever  surface  waves  (Lighthill,  1 9^2; 
Phillips,  19'5S)  that  provided  U in  this  frame  of  reference  docs 
not  vanish  at  some  level  (i.e.  provided  the  mean  flew  in  this 
frame  does  not  reverse,  in  which  case  we  have  a critical  layer) 
the  magnitude  of  these  corms  is  negligible  compared  with,  say, 
the  first  term  of  [2.6],  In  this  circumstance,  the  streamlines 
of  the  mean  flow  arc  the  same  as  they  would  be  if  the  flow  were 
non-turbulent  (though  not,  of  course,  irrotational  because  of 
the  mean  shear)  and  the  variations  in  surface  velocity  could  be 
calculated  ignoring  the  effects  of  the  •’■turbulence. 


Some  experimental  results  shed  light  on  this  question. 

The  same  Equation  f2.6]  is  relevant  to  the  air  flow  over  sur- 
face waves  where  also  there  is  shear,  turbulence  and  an  imposed 
wave  perturbation.  Laboratory  measurements  of  the  rate  of  de- 
crease with  height  of  the  wave  induced  fluctuations  have  been 
made  by  Karaki  and  Hsu  (ISao).  In  these  measurements,  there 
is  a critical  or  matched  layer  (where  the  wind  speed  U equals 
the  wave  speed  c)  quite  close  to  the  water  surface  so  that  the 
amplitude  and  phase  relations  are  distorted  between  the  sur- 
face profile  and  the  induced  air  flow  disturbance.  In  the 
region  above  the  critical  layer  however,  Karaki  and  Hsu  found 
that  the  (U,VJ)  perturbations  are,  within  experimental  accuracy, 
in  quadrature  (as  they  are  in  an  inviscid,  non-turbulent  flow). 
Furthermore,  the  decrease  of  the  amplitude  of  the  perturbations 
with  height  is  reasonably  consistent  (Figure  ) with  the 
Lighthill -Miles  approximation  for  the  perturbation  stream  function 
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V'  - ,(U(z) . * c)  c’kz  [2.8] 

which  is  valid  av/ay  from  the  critical  layer  in  a non-turbulcnt 
flovr. 

These  results  surest  that  insofar  as  the  wave-induced 
fluctuations  are  concerned,  the  turbulence  maintained  and 
generated  by  the  mean  flow  has  no  substantial  influence,  at 
least  in  the  air  flow  over  'waves.  Is  the  same  true  for  a cur- 
rent distribution  in  the  upper  ocean  moving  relative  to  an  in- 
ternal v:ave  field  supported  by  the  thermocline  below?  The 
principal  dynamical  difference  between  the  two  cases  concerns 
the  ratio  of  tv;o  time  scales  — the  readjustment  time  T of  the 

&L 

turbulence  to  a perturbation,  relative  to  the  time  that  a parcel 
of  fluid  taker,  to  move  from  one  wave  crest  to  the  next, 

X/[  U(s ) - c]  'Where  X is  the  wave  length.  If  this  ratio  is  large 
in  both  situations,  then  the  turbulence  has  little  opportunity 
to  adjust  to  a local  disturbance  and  its  structure  is  modified 
only  slightly.  If  it  is  small,  then  the  turbulence  is  approach- 
ing equilibrium  'with  the  locally  disturbed  mean  field  at  all 
stages . - - 

The  readjustment  time  of  the  turbulence  can,  very  plausibly, 
be  identified  'with  the  integral  time  scale  of  the  energy-con- 
taining eddies  in  a frame  cf  reference  moving  with  the  mean  flew  — 
this  is  the  appropriate  measure  of  the  "memory"  of  the  turbulent 
energetic  (and  Reynolds  stress)  structure.  There  is  now  good 
evidence  (Davies,  Fisher  and  Darratt,  19;>3)  that  in  shear  flow 
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this  time  scale  is  proportional  to  and  of  the  order  of  (SU/c^z)  ’ 
the  inverse  of  the  local  mean  shear,  so  that  the  ratio  of  the 
turbulent  adjustment  time  to  the  time  of  passage  over  one  wave 
length  is 


_ U(z)  - C 

Mdu/dz) 


[2.9] 


Tills  ratio,  of  course,  varies  in  each  flow  but  representative 
values  are  of  interest.  In  air  flow  over  waves  in  the  laboratory, 

such  as  in  Karahi  and  Hsu's  experiments,  above  the  critical  layer 

— ■ — -i 

U(z)  - C ~ 5 tn/sec,  X ~ 0.5  m and  dU/dz  ~ 50  sec  , decreasing 

with  height.  With  there  figures,  y ~ 0 .2  (increasing  with 
height)  so  that  in  this  region  the  turbulence  is  adjusting 
reasonably  rapidly  to  the  disturbed  flow.  One  would  therefore 
expect  lame r variations  in  Reynolds  stress  variations  <uw>, 

<u2>,  etc.,  than  when  ,•  » 1,  but  in  spite  of  this,  the  ’wave- 
induced  perturbations  seem  not  to  be  influenced  substantially 
(c.f.  Figure  ) . In  the  ocean,  with  a velocity  difference  U - C 
between  the  internal  'waves  and  the  water  near  the  surface  of 
0.5  m/sec,  an  internal  wave  wave: length  of  200  m and  a mean  eur- 


-2-1 

rent  shear  of  10  sec  , y *=  0.251  Representative  oceanic  values 


of  v are  apparently  of  the  same  order  as  laboratory  values;  we 
conclude  tentatively  that  the  perturbations  to  the  current  In- 
duced by  the  internal  waves  can  be  calculated  neglecting  the  in- 
fluence of  the  turbulence  to  an  accuracy  exemplified  by  the  com- 
parison of  Figure 
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2.3  Summary 

The  ocean  above  the  thernorlir.e  wili  have  mean  motions  on 
at  least  tv:o  distinct  scales: 

1.  Very  near  the  surface,  in  the  top  few  millimeters. 


there  is  a wind  drift  of  about-  3-df  percent  of  the  v/ind  speed  — 
laboratory  and  field  measurements  are  consistent  with  theoretical 
expectations.  It  is  disrupted  locally  by  wave  breaking,  but  is 
re-established  rapidly  by  the  v/ind.  It  responds  less  rapidly 
to  small  changes  in  wind  stress  and  will  tend  to  be  augmented 
somewhat  near  wave  crests  by  the  additional  wind  stress  (though, 
as  we  will  S'.e  later,  other  dynamical  wave -interact ion  effects 
are  usually  dominant). 


2.  Cn  a larger  depth  scale,  of  the  order  of  tens  of 


meters,  the  h or i rent 

al 

current  is  lik 

ely  to 

be 

variable  with 

depth  and  irregular; 

. j 

distortions 

of  this 

r» 

urront  distribu- 

tion produced  by  an 

inf 

ernal  wave  do 

der-end 

on 

tive  current  dir. 

tribation  itself  but  are  not  strongly  dependent  on  the  ambient 
turbulence  level. 
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3.  THE  VARIATIONS  IN  SURFACE  VELOCITY  RESULTING 
FROM  INTERNAL  WAVES 

We  are  interested  in  the  variations  in -the  velocity  dis- 

i 

tribution  near  the  ocean  surface  that  are  associated  with  in- 
ternal waves.  However,  the  small  scali  configuration  and 
properties  of  the  surface  are  coupled  to  longer  waves  that  are 

1 • i ■ 

• influenced  by  the -variations  in  the  phased  averaged  velocity 

i • . 

to  depths  of,  say,  10  m.  Consequently]  we  need  to  estimate  the 

* i 

variations  in  horizontal  current  near  the  surface  (in  this 

Sense)  that  are  induced  by  internal  waves,  but  we  are  not  here 
- ! 
concerned  with  the  variations  at  greater  depths. 

3-1  The  Irrotatior.nl  Solution 

When  there  is  no  mean  current,  or  the  current  is  uniform 
with  depth,  the  solution  is  elementary  and  well  known.  If 
d is  the  depth  #to  the  top  of  the  thermocline  (the  bottom  of 
the  mixed  layer)  and  a the  amplitude  of  the  wave  disturbance 
there,  then  in  a sinusoidal  wave  train  (or  a group  with  more 
than  two  or  tMjroe  waves)  the  velocity  induced  at  the  surface  is 
in  the  direction  of  wave  propagation  and  is  of  the  form 
(Phillips,  1963,  p.  167) 


when  the  vortical  displacement  at  depth  d is  C ^ a cos  (kx-nt). 

The  internal  wave  speed  C deDends  on  the  wave  length,  the 

| 

structure  of  the  thermocline  and  the  mode  number  of  the  inter- 
nal wave,  but  [3.1)  remains  valid  for  the  velocity  in  the  mixed 


Wp 

t 
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layer.  For  a sharp  thermocline  in  deep  water,  for  the  lowest 
mode , 


c , H « / Slfi.  1 

k ) p k(l  + coth  kd) 


J 


[3-2] 


For  internal  waves  that  are  long  compared  with  the  thermocline 
depth  d,  kd  « 1 and 


l. 


[3-31 


Characteristically,  for  a thermocline  depth  of  50  m and  a 

-3 

density  jump  6 p/p  ~ 10  , C ^ 0.7  m/sec  for  internal  wave 
lengths  of  about  200  m or  more ; smaller  for  shorter  wave  lengths. 
The  magnitude  of  the  orbital  velocity  differences  induced  by  the 
internal  wave  near  the  surface  is  2an;  if  the  wave  amplitude  is 
5 m then  (with  the  figures  above)  this  difference  is  only  0.2 
m/sec  for  a 200  m internal  wave,  decreasing  in  inverse  propor- 
tion as  the  internal  wave  length  A increases.  The  maximum 

-3  -1 

strain  rate  near  the  surface  is  ank  = 3 >-'  10  sec  for  a 200  m 

• -2 

wave;  this  is  proportional  to  A 

Fluid  velocities  associated  with  surface  wind-generated 
waves  are  usually  larger  than  the  value  calculated  above  and 
(on  account  of  the  smaller  scales  of  surface  waves)  the  strain 
rates  are  very  much  larger.  Clearly,  the  motions  induced  by  the 
internal  waves  are  generally  a very  small  perturbation  (but  of 
relatively  large  scale)  on  an  active  and  noisy  background;  the 


r-‘ 
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possibility  of  detecting  the  signal . depends  on  (a)  whether  the 
velocity  variations  induced  near  the  surface  by  the  internal 
wave  can*  under  certain  natural  conditions,  be  augmented  and 
(b)  the  existence  of  one  or  more  indicators  (or  dynamical  pro- 
cesses) that  are  sufficiently  responsive  to  these  variations. 

3.2  Interactions  with  Surface  Currents 

The  express  ion  [ 3.1]  holds  only  - in  the  absence  of  any  ver- 
tical variation  in  the  m-aan  horizontal'  current.  When  there  is 
structure  in  the  current,  the  surface  velocity  variations  pro- 
duced by  internal  waves  can  be  modified  substantially  as  the 
following  analysis  demonstrates . 

In  Section  2 of  this  work,  it  was  shown  that  the  wave  in- 
duced perturbations  appear  to  be  substantially  unaffected  by 
the  turbulence  provided  there  is  no  critical  layer  where  the 
component  of  the  mean  current  velocity,  in  the  direction  of 
propagation  of  the  internal  wave  is  equal  to  the  internal  wave 
speed.  Under  these  circumstances,  [2.G]  reduces  to 


u a + v f 

dx  dZ 


[3< 


~for~the~~c ornpoheht  " of  'phase "’averaged  vortxcity  (the  mean  plus 
the  wave-induced  variation)  normal  to  the  direction  cf  wave 
propagation.  Since  fi  - SU/hz  - 9W/Bx,  this  equation  is 
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or,  after  a little  manipulation  and  use  of  the  incompressibility 
condition  dU/dx  + SW/dz  = 0 (note  b/by  = 0), 


JL  TT  £U  ^ JL  / d u 

az  u ax  * az  \ A | az 


j_  (,.  sw  i 
ax  J ax  I 


= 0 


[3.5] 


Integrate,  this  equation  vertically  from  the  top  of  the  disturbed 
thermocline  at  z = -d(x)  to  the  surface: 


.A 

ax 


w 


au  aw), 
az  “ ax*id 
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[3-6] 


and  compare  the  order  of  magnitude  of  these  terms. 

In  the  frame  of  reference  moving  with  the  waves,  the  mean 
velocity  U Is  of  the  order  C in  magnitude  and  the  variations  in 
it  are  of  order  Ca/A.  By  the  incompressibility  condition,  for 
internal  waves  appreciably  longer  than  the  thermocline  depth  d, 

W ~ C (a/A)(d/A).  Now  the  third  term  contains  the  quantity 
dU/dz  which  is  either  of  order  Ca/Ad  or  dominated  by  the  ambient 
mean  shear  at  the  depth  -d.  In  order  to  refer  the  mean  surface 
current  to  the  current  speed  (if  any)  at  the  thermocline,  we 
will  need  to  choose  d at  a level  where  the  mean  shear  is  small 
so  that  the  difference  is  clearly  defined;  the  second  alterna- 
tive is  therefore  excluded.  Consequently  the  orders  of  magni- 
tude are: 
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Relative  Order 

1st  two  terms 

C2a/A2 

1 

Third  term 

C2 (a/A )2 (d/A) 2 . l/A 

(a/A) (d/A)2 

Fourth  term 

C2 (a/A) (d/A)2  1/A 

(d/A)2 

Now  a/A  « 1 (generally  less  than  0.2)  and  in  many  cases  of 
interest  d/A  <0.3  approximately.  In  the  long  wave  length  limit 
as  (d/A)2  -*-0,  Equation  [3*6]  reduces  simply  to 


or 


s 

U 2 - U . * const  w.r.t.  x. 

o -d 


[3*7] 


It  can  also  be  shewn,  by  taking  the  component  of  the  vorticity 
equation  in  the  direction  of  wave  propagation  that  the  surface 
current  transverse  to  the  direction  of  wave  propagation  is  un- 
affected by  the  internal  'waves : 


v **  const, 
o 


[3.81 


The  detailed  calculation  can  be  extracted  from  Phillips  and 
3anner  (1971!). 

Note  that  the  expression  [3*7]  is  valid  only  in  the  long 
wave  length  domain,  when  (A/d)2  » 1.  The  case  of  shorter 
wave  lengths  can  be  considered  numerically  by  integration  of 
Equation  [3*1-]  out  this  work  is  not  yet  complete.  Note  also 
that  [3*7],  relating  the  variations  in  surface  current  to  the 
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velocities  at  the  thermocline  depth,  does  not  Involve  the  de- 
tails of  the  current  structure  in  between.  It  is  an  inte- 
grated form  of  the  vortioity  conservation  equation;  to  this 
order  in  (Vd)  the  vorticity  is  simply  the  mean  shear  whose 
integral  is  the  velocity  difference  across  the  layer. 

3.3  Surface  Current  Modulations 

1 

i 

The  velocities  U and  U , are  taken  in  a frame  of  reference 
o -d 

moving  with  the  internal  wave,  so  that  if  u(x)  represents  the 

i ' 

orbital  velocity  at  the  thermocline  and  Ue (x)  the  surface  cur- 
rent (relative  to  the  datum  against  which  C is  measured),  then 
[3-7]  becomes 

[U  (x)  - CJ2  - [ u (x ) - C]2  * const  [3-9] 

s 

Now  in  a periodic  internal  wave,  u(x)  = 0 at  the  phase  of  the 
wave  'where  the  thermocline  displacement  vanishes  and  if  Ug  is 
the  surface  current  at  this  phase  point,  then 

[Us(x)  - C]2  - (u(x)  - C]2  * [Us  - C]2  - C2  • 

and 

Us(x)  = C - [ (C-u(x)  )2  - Ug  (2C  - Ug)]2  [3.10] 

the  negative  sign  being  chosen  since  U - C < 0 when  there  is  no 

critical  layer.  Note  that  when  the  surface  current  vanishes, 

U = 0 and  [3-10]  reduces  to  U (x)  = u(x),  the  long  wave  length 
s s 

limit  for  internal  waves  in  an  irrotational  flow.  Since 


HYDRONAUTICS,  Incorporated 

-24- 

u = - (C  a/d)  cos  x for  long  Internal  waves,  where  x Is  the 
phase  of  the  wave,  x - 0 corresponding  to  an  internal  wave 
crest,  [3*10]  can  be  written 

U8(x)  . 1 

~~~  = 1 - [(1  + (a/d)  cos  x)£  - r(2-r)p  [3*11] 

where  r = U /C.  If  U and  therefore  r are  positive,  i.e.  when 

S S 

the  mean  surface  current  has  a component  in  the  same  direction 
as  that  of  the  internal  wave  propagation,  the  modulations  in 
surface  current  are  greater  than  in  a purely  irrotational  wave. 
Us(x)  is  greatest  (in  the  positive  x-direction)  when  cos  x = -1* 
above  the  internal  wave  troughs. 

Figure  5 shows  some  representative  distributions  of  surface 
current  with  respect  to  phase  of  the  internal  wave  for  various 
values  of  r when  (a/d)  = 0.1.  Note  that  the  amplitude  of  the 
modulations  in  surface  current  increases  for  a fixed  a/d  as  the 
mean  surface  current  increases;  when  r *=  u /c  = 0.5641  (for 
a/d  = 0.1)  the  magnitude  of  the  surface  current  above  the  in- 
ternal wave  trough  becomes  equal  to  the  phase  speed  of  the 
internal  wave.  For  larger  values  of  U /C  the  water  here  is  in 
fact  surging  ahead  faster  than  the  wave  and,  in  a frame  of 
reference  moving  with  tne  wave,  there  is  a closed  eddy  in  the 
mean  flow  near  the  surface.  The  detailed  calculation  becomes 
unreliable  under  these  circumstances  (Reynolds  stress  varia- 
tions being  important)  but  the  pattern  of  mean  streamlines  be- 

co:n>  s qualitatively  as  in  Figure  6.  When  U /C  is  substantially 

s 
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greater  than  1,  there  is  a submerged  critical  layer,  with 
energy  exchange  between  the  current  and  the  internal  waves. 

The  magnitude  of  the  modulations  in  surface  current  as 
functions  of  r and  a/d  is  shown  in  Figure  7 ; the  amplifying 
effect  of  a surface  current  moving  with  the  internal  waves  is 
evident.  The  curves  for  various  values  of  a/d  terminate  at 
the  point  where  the  velocity  of  the  water  at  the  surface  above 
the  Internal  wave  trough  is  equal  to  the  phase  speed  of  the 
wave.  As  Figure  5 shows,  the  forward  surging  motion  is  both 
more  localized  and  mere  intense  than  the  backwards  surge  over 
the  internal  wave  crests. 

Similar  information  is  given  in  Figure  8.  Here,  for 
example,  an  internal  wave  with  amplitude  a = 0.075  d in  a sur- 
face current  with  mean  speed  0 . -1 8 C has  surface  current  modu- 
lations twice  that  of  an  irrctational  wave  with  the  same  amp- 
litude and  speed,  the  velocity  difference  at  the  surface  being 
0.3  C between  points  over  the  internal  wave  trough  and  crest. 
Points  that  lie  in  the  upper  right,  outside  the  contours, 
identify  conditions  in  which  a closed  eddy  exists  at  the  sur- 
face vrith  . fluid  over  the  trough  surging  ahead  faster  than  the 
wave  profile. 

When  the  surface  current  is  inclined  to  the  direction  of 

propagation  of  the  internal  wave,  the  component  of  the  current 

normal  to  the  propagation  direction  is  not  affected  (Phillips 

and  Banner,  157- ).  In  a frame  of  reference  moving  with  the 

waves,  the  pattern  cf  surface  streamlines  is  as  shown  in  the 

upper  part  of  Figure  9 for  the  case  (U  ) /C  =0.5,  a/d  = 0.1 

s x 
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with  an  angle  of  -*f 3°  between  the  surface  current  and  the 
propagation  direction  of  the  waves*  When  (U  ) /C  becomes 
larger. (or  the  wave  amplitude  a/d  increases)  the  flow  reverses 
in  the  detached  eddy  end  as  indicated  qualitatively  in  Fig- 
ure 9b,  a divergence  line  and  a convergence  lino  appears.  The 
latter  (which  will  be  seen  later  is  significant  for  slick 
formation)  lies  between  an  internal  wave  crest  and  the  fol- 
lowing trough.  As  the  x -component  of  U increases,  the  two 
lines  separate  towards  adjacent  crests  until,  at  a greater 
value  of  U /C  > 1,  they  coalesce  at  the  crest  and  disappear  — 
at  all  phases  of  the  internal  wave , the  surface  current  is  then 
overtaking  the  wave.  Figure  10  illustrates  the  various  domains 
of  behavior;  the  right-hand  branch  being  rather  unreliable 
since  a critical  layer  is  then  developing  and  the  essentially 
non-turoulent  calculation  here  becomes  unreliable. 

3 Ji  Summary 

1.  The  modulations  in  surface  curreiit  produced  by 
internal  waves  are  amplified  when  a mean  surface  current  is 
moving  in  the  direction  of  internal  wave  tronamnt Ion,  suppressed 
when  moving  against. 

2.  The  amplification  ratio  becomes  large  'when  the 
component  of  mean  current  speed  along  C,  Ug } is  close  to  the 
propagation  speed  C of  the  internal  waves  and  the  internal  wave 
slope  is  small. 
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3.  When  U /C  is  near  unity,  the  surface  current  re- 
verses  direction  relative  to  the  internal  wave  pattern,  even 
for  internal  waves  of  small  slope.  This  forms  convergence  lines 
between  a crest  and  the  following  trough  for  slicks  and  other 
surface  material.  Above  the  internal  wave  trough,  the  water 
near  the  surface  moves  faster, than  the  internal  wave  pattern. 

•’(.  Thesp  effects  are  strongest  when  the  internal 
wave  length  is  much  greater  than  the  thcrmocline  depth.  Cal- 
culations for  cases  when  the  two  are  comparable  are  presently 
being  done.  • 
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4.  SURFACE  FILMS 

4.1  General  Properties 

Surface  films  are  very  widespread  on  the  ocean  surface, 
whether  as  dispersed  or  saturated  monolayers  or  thicker  films, 
and  seem  generally  to  be  of  biological  origin.  Anyone  who  has 
worked  with  precise  wave  experiments  in  the  laboratory  can  at- 
test to  the  difficulty  of  getting  rid  of  surface  films,  even 
on  ’’clean"  chlorinated  and  filtered  water.  A general  account 
of  the  phenomena  of  surface  films  is  given  by  Davies  and  Rideal 

(19o3).  The  thickness  of  a monolayer  is  of  the  order  10  craj 

-3  -1 

it  has  a surface  viscosity  of  10  J to  1 e.g.s.  units  (gm  sec  ) 

and  its  compressibility  is  very  non-linear  with  film  density. 

Conservation  equations  can  be  written  for  the  surface  film 
as  follows.  If  v is  the  film  velocity  over  the  surface,  y the 
film  density  (mass/area),  then  by  mass  conservation 


+ V • (yv)  - 0 


[4.1] 


+ y y-v  = 0 


[4.2] 


where  the  divergence  is,  of  course,  on  the  surface.  The  momentum 
equation  for  the  film  is 


(y  v)  * 


+ u.ysv  + F 


[4.3] 
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where  p^,  = p^y)  Is  the  film  pressure,  d the  surface  viscosity 
and  P the  drag  force  per  unit  area  resulting  from  (a)  the  wind 
stress  at  the  top  and  (b)  the  lower  stress  from  any  relative  motion 

of  the  film  in  the  underlying  fluid.  Since  F may  be  of  the 

-7 

order  1 dyne/cm2  while  y ~ 10  gm  cm2  the  acceleration  terms 
in  [t.3]  are  negligible  and  the  stress  balance  is,  in  essence, 
static: 

dpf 

. * xjy  + 4V2V  + F J=  0 [4.4] 

The  quantity  Spf/oY  is  related  to  the  two-dimensional  compressi- 
bility modulus  c^  ^ (Davies  and  Rideal,  p.  263) 

cs*1  “ y Spf/cV  [4‘5] 

For  a very  disperse  film,  c — pr.  -*  0 and  u - 0 as  y -*•  0. 

S -1A  2 -1 

For  tightly  packed  monolayers,  c ~ 10  dynes  cm  for  liquid 
-11  S 

films  and  c ~ 10  for  fatty  acids.  Other  figures  and  ex- 
s 

perimental  results  are  given  in  the  book  cited. 

In  the  ocean,  dispersed  films  are  not  usually  visible  to 
the  eye  while  concentrated  films  can  be  seen  generally  by  their 
influence  on  short  wave  breaking  (as  will  be  discussed  later). 

Films  shewing  interference  coloring  are  much  thicker  than 
monolayers  —generally  oil  spills. 
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4,2  Influence  of  Short  Waves  on  Average  Film  Density 

The  most  vigorous  motions  at  the  ocean  'surface  result  from 
surface  waves,  but  the  length  scale  of  this  motion  is  small 
compared  with  internal  wave  lengths  of  interest.  Accordingly, 
let  i£(x,t)  represent  the  film  velocity  averaged  over  several 
wave  lengths  of  the  short  surface  wave  field,  and  v'(x,t)  repre- 
sent the  fluctuations  about  this  average  produced  by  small  scale 
surface  waves.  Similarly,  let  T(x,t)  and  y’(x,t)  represent  the 
mean  and  fluctuating  parts  of  the  film  density. 

Equation  [4.1],  when  locally  averaged  in  this  way,  becomes 

+ V*  (TV  + yTyf)  - 0 [4.6] 

Since  the  root  mean  square  wave-induced  velocity  (v2)?  is 
generally  large  compared  with  V,  the  influence  of  the  second- 
term  on  the  mean  film  density  must  be  estimated. 

Consider  a wave  moving  across  the  surface  where  there  is 
a diffuse  film  of  density  yQ.  If  the  wave  profile  is 

C = a cos  (k*x  - nt) 


then  the  surface  velocity  in  deep  water  is 


~ y*r  » ' r lyr1'  ** 
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where  1 is  a unit  vector  in  the  direction  of  k and  the  fre- 
.1 

quency  n = (gk)‘?.  In  a frame  of  reference  moving  with  speed 
x 

c = of  the  wave 


v * an  coskx  - c 


C A.8] 


where  the  x-direction  is  taken  in  the  direction  of  v/ave  travel. 
The  conservation  Equation  [4,1]  becomes 


^'/[an  coskx  * c]  *)  « 0 

since  in  this  frame  of  reference,  the  motion  is  steady.  Thus 


y [an  coskx  - c]  « const  = - F^c  [4-9] 


Nov;  for  surface  waves  c » an,  so  that 


-1 


y « r ( 1 - — coskx  \ a f / 1 + — coskx  + 0 [ — 

O i C ( O \ C C j 


*\ 


and  the  convariance  between  the  fluctuations  y’  and  v* 


y ' v ' *=  a2na/e 


[4.10] 


With  a spectrum  of  short  wques"  ,*{n),  essentially  uni-directional 


y1  v’ 


6 


n3'i  (n)  dn 


P ’ 7 c*  *■>  ■*  r 

o w * ^ 


[4.11] 


[4.12] 
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where  Vw  is  the  effective  surface  film  flux  velocity  induced 
by  the  waves  — in  fact,  the  Stokes  drift.  In  a saturated 
spectrum, i(n)  * ^g2n~s,  n > nQ, 

Vw  • ^s/no  [ ^ .13] 

where  n^  (rad/sec)  is  the  frequency  of  the  spectral  peak  and 

w _p  | 

£ ^ 1.1  x 10  . The  drift  velocity  V is  thus  approximately 

W I | 

1 percent  of  the  speed  of  the  highest  waves . The* mean  local 
film  density  I’  is  therefore  governed  by 

f£  * V-  (r(V  + . 0 t4.»] 

4 . 3 Variations  in  Average  Film  Density  Produced  by  internal  - 
Waves 

In  a dispersed  film,  c and  u are  both  small  and  [4.4]  t 

s 

reduces  to  F --  C;  the  local  wind  stress  at  the  upper  surface 

of  the  film  balances  the  drag  of  the  water  below.  The  film 

simply  moves  with  the  free  surface  at  a velocity  that  is  the  'i 
i v > 

vector  sum  of  the  'surface’  current  U„  discussed  in  Section  3> 

■TV'S 

the  Stokes  drift  V and  microscale  surface  wind  drift  q the 

'■c 

last  of  which  is  in  the  direction  cf  the  wind  stress  with 
magnitude  approximately  3 or  4 percent  of  the  wind  speed  U . 

The  conservation  equation  for  the  mean  film  density  [4.14]  is 
therefore 

f + 7-  rr(ys  + + 9)1  * o 


[4.15] 


£ 
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Take  coordinates  with  x in  the  direction  of  wave  propagation 
and  moving  with  the  speed  C of  the  internal  waves,  in  which 
frame  [ -4.15]  becomes 

7 f + Ir  [r(us  + ’ + vw  - c>x]  = 0 t».:6] 

involving  the  components  of  the  velocities  in  the  x-direction. 

In  a very  long  internal  wave  train  a steady  state  may  be  set 
up  in  which 


I?  !r<Us  + ■>  + v„  - °>x]  * 0 

so  that  the  average  local  film  density  varies  inversely  with 
the  x-component  of  the  surface  velocity  including  the  wind  drift 
observed  as  we  move  with  the  internal  wave.  If  at  some  phase 
of  the  wave 

Us(x)  + q + Vw  - C - 0 [4.18] 


the  solution  fails ; f -*•  » and  the  film  is  clearly  no  longer 
diffuse. 

The  development  of  patches  of  high  film  concentration 
when  a pulse  of  internal  waves  enters  a region  of  diffuse  film 
can  be  shown  simply.  From  [4.16] 


F § + <us  - c + i + V 


i ar 
r ax 


5Us(x) 
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since  C and  q are  constant;  we  neglect  the  variation  in  V on 

W 

the  ground  that  it  is  dominated  by  the  longest  waves  (c.f.  4.13). 
Thus 


d r 

Dt  in  r 


[4.20] 


where  D/Dt  represents  the  derivative  following  the  motion  of 

the  film  and  V is  the  initial  film  density.  If  at  some  phase 

o 

of  the  wave  [4.17]  is  satisfied,  a convergence  line  exists  in 
which  5U  /?;•;  < 0;  material  points  of  the  film  will  move  towards 
this  line  and  (from  4,20)  the  film  density  increases  exponentially 


r(4 ) — exp  / s (t)dt 


where  s(t)  is  the  convergence  - oU/Sx  experienced  by  the  material 
point.  As  t increases,  s(t)  -*■  s(0),  the  convergence  at_  the  con- 
vergence line  and 


r(4 ) —-ro  exp  [ s (0) t] . 


[4.21] 


The  film  density  therefore  increases  exponentially  with  time 
until  it  saturates > even  if  the  initial  film  is  very -disperse  — 
the  lower  the  initial  film  density  the  longer  the  film  will  take 
to  saturate  and  the  narrower  the  ultimate  saturated  region,  ex- 
cept at  the  leading  edge . of  an  internal  wave  group  where  material 
will  continue  to  be  collected. 
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The  development  of  concentrated  slicks  is  therefore  de- 
pendent on  [ 4 . 18]  being  satisfied  at  some  phase  of  the  internal 
wave.  The  condition  is  illustrated  in  Figure  11.  If  Vw  + q * 0 
(say  if  the  wind  is  blowing  along  the  crests  of  the  internal 

waves)  slicks  are  formed  for  values  of  a/d  and  U /C  that  cor- 

s 

respond  to  points  inside  the  bowl -shaped  curve.  If  there  is  no 
mean  surface  current  component  in  the  direction  of  wave  travel, 
the  amplification  effect  on  the  surface  speed  modulations  dis- 
appears and  slicks  are  formed  only  inside  the  V-shaped  region 
shown.  In  any  event,  for  slicks  to  be  associated  with  internal 
waves,  it  is  necessary  that  the  vector  sum  of  the  average  cur- 
rent, the  wind  drift  (about  3 percent  of  the  wind  velocity)  and 
the  wave  drift  (about  1 percent  of  the  speed  of  the  highest  waves) 
have  a component  in  the  direction  of  internal  wave  travel  that 
is  of  the  same  order  as  (though  not  necessarily  close  to)  the 
speed  of  the  internal  waves.  The  range  of  conditions  under  which 
they  will  be  seen  increases  with  the  ratio  of  internal  wave  am- 
plitude to  thermocline  depth. 

If  the  average  density  of  surface  film  is  very  low,  the 
material  gatb®»v?d  at  the  convergence  lines  will  lie  in  a thin 
strip  and  may  be  dispersed  again  by  breaking  waves  incident  on 
this  region.  The  influence  of  saturated  or  nearly  saturated 
slicks  on  short  wave  breaking  will  be  considered  in  Section  6 of 
this  report.  Their  influence  on  the  damping  of  unbroken  waves 
is .well  known  (see,  for  example,  Phillips,  1968,  p.  37).  The 
suppression  of  short  wave  breaking  together  with  the  increased 
attenuation-  of  capillary  waves  gives  the  surface  of  an  extended 


-36- 


HYDRONAUTICS,  Incorporated 


saturated  slick  a glassy  (or  oilyl)  appearance.  LaFond’s  (1962) 
early  observations  on  the  associations  between  internal  waves 
and  moving  slicks  were  concerned  with  extended  saturated  slicks, 
easily  visible  by  eye.  The  associations  between  the  two  phe- 
nomena are  clear  in  h.^s  observations,  though  his  simple  explana- 
tion is  untenable. 

If  the  condition  [4.18]  is  not  satisfied,  that  is,  if  the 
oceanic  conditions  are  represented  by  points  outside  the  domains 
of  Figure  11,  variations  in  film  density  will  be  generated  by 
internal  waves,  but  whether  or  not  a slick  is  formed  depends 
upon  the  average  ambient  density  of  the  surface  film.  For, 
from  [ . 17  ] the  local  mean  film  density  T(x)  is  given  by 

r(x)  [Ux(x)  + q + Vw  - C]  = r (Us  + q + Vw  - C]  [4.22] 

IfU  +q+V  - C > 0,  the  maximum  film  density  will  occur 
s w 

when  U is  least,  above  the  internal  wave  erest  (see  Figure  5)* 
s 

On  the  other  hand,  ifU  + q + V - C < 0,  the  internal  wave 

s w 

is  overtaking  the  surface  material  everywhere  and  the  maximum 
film  density  will  occur  where  U is  greatest,  above  an  internal 

O 

wave  trough.  We  see  later  that  if  the  local  mean  surface  film 
is  unsaturated  (not  tightly  packed,  but  fairly  close  to  it) 
variations  in  mean  film  density  can  have  a significant  effect 
on  wave  breaking. 
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4.4  Summary 

1.  Internal  waves  modulate  the  density  of  naturally 

occurring  surface  films;  this  in  turn  can  Influence  the  break- 

j ■ 

ing  of  small-scale  waves  (see  Section  6.2)  even  when  the  average 
film  density  is  such  that  no  slick, is  visible. 


2.  Visible  slicks  are  formed  when  the  vector  sum  of 
the  surface  current,  the  wind  drifil  i (about  3 percent  of  the  wind 
velocity)  and  the  wave  drift  about  1 percent  of  the  speed  of 
the  highest  waves ) has  a component  in  the  .direction  of  wave1 
travel  that  is  of  the  same  order  as  (though  not  necessarily 
very  close  to)  the  speed  of  the  internal  wave.  The  range  of 
conditions  under  which  visible  slicks  will  be  seen  increases 
with  the  ratio  of  internal  wave  amplitude  to  thermocline  depth. 
The  convergence  line  associated  with  the  slick  lies  between  an 

■ P 

internal  wave  crest  and  the  following  through,  if  U + q + V < C; 

s w 

between  a trough  and  the  following  crest  if  U + q - V.  > C. 
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5.  ENERGY  MODULATIONS  IN  A SURFACE  WAVE  TRAIN 
PRODUCED  BY  DIRECT  INTERACTION' 

The  direct  Interaction  between  surface  waves  and  much  longer 
internal  waves  produces  modulations  in  energy  density  of  the  sur- 
face waves  that  can  certainly  be  observed  under  ideal  conditions 
in  the  laboratory.  Whether  or  not  modulations  are  generated  in 
the  field  in  this  way  depends  on  whether  other  effects  with  a 
shorter  time  scale  than  the  surface-internal  wave  interaction 
will  not,  or  will  overwhelm  this  process.  T^ere  are  two  ways  of 
considering  the  dynamics  of  the  simple  direct  interaction;  they 
are  complimentary  and  each  useful  in  appropriate  circumstances. 

It  can  be  considered  from  the  point  of  view  of  a wave  train  on  a 
slowly  varying  current  or  from  the  viewpoint  of  resonant  inter- 
actions among  internal  and  surface  waves. 

5*1  Waves  on  a Slowly  Varying  Current 

The  energy  density  E of  a wave  train  in  a slowly  varying 
current  U(x,t)  is  specified  by  the  equation 

) > +•  S • 7U  * 0 * [5*1] 

where  c is  the  local  grout)  velocity  and  S the  radiation  stress 
tensor  (integrated  Reynolds  stress),  which  in  deep  water  takes 
the  form 


+ 


c 


[5.2] 


where  the  1-dircction  is  that  of  the  local  wave  propagation.  The 
pattern  of  wave-numbers  is  given  by  the  kinematic  equation 
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ak 

7 (0  + l *£)  - O [5.3] 


where  o = c(k)  is  the  intrinsic  frequency;  for  surface  waves 
a - (gk)^.  The  derivation  of  these  equations  is  given  by 
Phillips  (1968). 

When  a pattern  of  internal  waves  (either  a continuous  train 
or  a steady  pulse)  is  propagating  v/ith  speed  C beneath  the  sur- 
face, a distribution  of  surface  current  of  the  form  Ug(x-Ct)  is 
set  up  (with  mean  U ) and  this  interacts  with  the  surface  waves. 
In  the  absence  of  effects  such  as  wind  input,  such  as  wave  break- 
ing or  attenuation  (which  will  be  considered  later)  the  distri- 
bution of  surface  wave  energy  can  be  found  by  integration  of 
[5.1]  and  [5 -3].  In  some  circumstances,  if  we  move  with  the 
internal  waves,  a steady  state  disturbance  of  surface  wave 

energy  is  established  and  solutions  for  U =0,  no  surface  cur- 

s 

rent,  have  been  given  by  Phillips  (1971)  and  Gargett  and  Hughes 

(1972).  These  results  are  generalized  below.  When  the  Internal 

and  surface  waves  are  collinear,  the  distribution  of  local  sur- 

1. 

face  wave  intrinsic  speed  c * (g/k)^  is  given  by  the  solution 
to  the  quadratic 


c*  c + Us(x)  - C 

c — _ 

0 c + U - C 
o s 


[5.*] 


where  cq  is  the  speed  at  the  point  where  Ug 


U . Consequently 
s 
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c(x) 

c ” 


c + (x)  - C)(c  + Ue  - C)}^ 


2(c  + U 

' o s 


C) 


The  alternative  root  is  eliminated  since  as  U u , c - 

s s 

The  distribution  of  energy  density  in  the  surface  wave  is 


c . 
o 


. c (u  - c + ie  ) 

eo  = c(x)  (us(x)  - c + foot 

Note  the  singularity  in  this  solution  if 

Ug(x)  + £c(x)  - C - 0 

when  the  energy  flux  velocity  U + (1/2 )c  of  the  surface  wave  is 

s 

just  equal  to  the  nhase  speed  of  the  internal  waves.  At  this 
point  the  surface  wave  speed  can  be  shown  (Fhillips,  1971)  to 
become  complex  and  waves  can  no  longer  propagate  past  this  point. 
This  is  the  phenomenon  of  blockage  — surface  wave  energy  tends 
to  accumulate  at  particular  phase  points  of  the  Internal  wave, 
and  the  way  that  it  accumulates  can  be  described  by  an  unsteady 
analysis  based  on  [5.1].  The  modulations  In  surface  velocity 
required  to  produce  blockage  are  illustrated  in  Figure  12.  Note 
that,  at  the  points  where  the  curve  touches  the  horizontal  axis, 
a very  small  perturbation  in  a uniform  current  will  prod\ice 
"blockage”  — this  is  the  phenomenon  of  resonance  with  Internal 
and  surface  waves  of  small  slope,  described  by  Lewis,  Lake  and 
Ko  (1974)  for  Us  = 0 (no  mean  surface  current).  The  condition 
for  resonance  is  that 


[5-7] 


[5-6] 
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[5.8] 

For  a surface  current  in  the  direction  of  travel  of  the  wave, 
resonance  occurs  at  a shorter  wave  length  than  in  the  absence  of 
a current^  if  it  is  in  the  opposite  direction,  resonant  surface 
waves  are  longer.  It  will  be  seen  later  that  at  short  wave 
lengths,  about  15  centimeters  or  less,  resonance  or  blockage  in 
the  presence  of  wind  under  natural  conditions  is  probably  not  a 
significant  dynamical  process  being  overwhelmed  by  other  effects 
such  as  wind  energy  input  and  wave  breaking. 

5.2  Non-Resonance  Modulations  in  Longer  Waves 

Longer  wave s or  swell  that  move  considerably  faster  than 
the  internal  wove s may  not  suffer  blockage  or  undergo  resonance, 
but  they  are  modulated  by  the  internal  wave,  particularly  when 
a surface  current  is  involved.  It  will  be  found  later  that  im- 
portant effects  can  result  from  these  modulations  by  further 
interaction  with  short  breaking  waves. 

Consider  first  the  case  when  the  surface  current,  surface 
wave  and  internal  wave  are  all  collinoar.  From  [5-4]  for  dif- 
ferential increments  5U  in  surface  current 


U + - C = 0 

s 2 o 


so  that 


25c  = 5c  + 5U 

c0  c + U - C 
o s 


c 

o 


c + 2U  - 
o s 


5U 


6c  « 


2C 


[5.9] 
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and  from  [5.6] 


Ec  [U  (x)  + £c(x)  - C]  * const 
s 


5E  be 
E c 


6U  + 4&c 
s z 

“s  + i'c  - c 


XU 


0 


With  the  aid  of  [5-9]  and  after  a little  algebra,  it  is  found 
that 


'For  long  surface  waves  travelling  much  faster  than  the  internal 
waves,  |y|  « 1 and  the  modulations  in  energy  density  are 
approximately 


[5.H] 


These  results  can  be  generalized  to  situations  in  which  the  in- 
ternal wave  motion,  the  surface  current  and  the  surface  wave  are 
not  all  collinear.  Consider  the  geometry  specified  by  the  fol- 
lowing diagram 


W|M 
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Taking  a reference  frame  moving  in  the  x-direction  with  the  in- 
ternal wave  (as  before)  the  kinematics  of  the  surface  wave  are 
specified  by  [5-3]  as 

•k  to  + i * 0 

or 

+ k[Ug(x)  cos  9 + Vs  sin  0 - C cos  8]}  » 0 [5.121 

and  from  v x k = 0,  we  have 

» (k  £in  9>  [5.13] 

Accord ingly } the  variations  in  wave  pattern  are  given  by 
c 5k  + 6k  C[Us(x)  - C]  cos  0 + Vs  sin  0}  + 6U  k cos  8 
“ k t(us(x)  * C)  sin  9 - Vp  cos  0]  60  = 0 


.... 
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6k  sin  fl  + S0  k cos  8*0 

Solving*  -the  variation  In  wave  number  magnitude  is  given  by 


SL  . . - £2UL 6U 

k +,  U - C us 

i s 


and  £he  variation  in  direction  by 

I | 

6S  * Vo  vC'  ?v  C ros  [5.15] 

; s 

i 

i 

involving  only  the  component  of  surface  current  In  the  x-direction. 
Since  Sc/cq  * - |(fk/k),  this  reduces  to  [5-9]  when  8 = 0. 

The  energy  density  distribution  is  given  by 


[Ec  (-C  + U + £c  cos  8)}  * 0 
oX  s 


[5.16] 


(This  generalization  of  [5.6]  Is  obtained  most  readily  from 
Bretherton  and  Garrett's  action  conservation  principle.)  Thus 
variations  are  given  by 

, 6U  + J(6c)  cos  0 - ic  sin  0 68 

~ + ££  + — 5 , = o 

E c ic  cos  8 - C + U 

s 

From  this*  and  with  the  aid  of  [5.14]  and  [5. 15]  it  is  found 
after  some  algebra  that 


6E 


* -2 


6U 


1 + cos  2 


c f (cos  8 


9 + y(2  + cos  0) 


+ 2y ) (1  + 2y) 


[5-17] 
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where  y * (U  - C/c).  When  0 = 0,  this  reduces  to  [5. 10].  For 
long  surface  wave  travelling  much  faster  than  either  the  internal 
waves  or  the  surface  current,  y « 1 and 


tE  ^ (1  4-  cos  5Us 

E cos  6 c 


6U 

-4  cos  S — — 
c 


[5.18] 


The  simplicity  of  this  generalization  of  [5. 11]  is  rather  sur- 
prising. These  results  will  be  used  later. 


5.3  Resonance 

Three  wave  trains  of  small  slope  can  undergo  resonant  in- 
teractions, exchanging  substantial  amounts  of  energy  amongst 
themselves  if  their  wave-numbers  k.  and  frequencies  n simul- 
taneously  obey  the  relations 


*1  “ £a  - k3 

nx  - n2  * n3 


[5.19] 


Suppose  the  v.'aVe  trains  1 and  2 represent  surface  waves  of 
almost  the  same  wavelength,  so  that  kx  * k2  + 6k,  and  let  the 
wave  train  3 represent  an  internal  wave . With  a surface  current 
U , assumed  constant,  the  frequencies  na  =*  cjs  + k2  * U , 
n,  = aa  + 6a  + (k3  + 6 k)  • U , so  that  the  conditions  [5.13] 
become 

6k  - k3  - . 


= n. 


5a  + 6k  • U 
~ ~s 


[5.20] 
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where,  as  before,  o represents  the  intrinsic  frequency.  From 
the  first  of  [5.20]  resonance  can  occur  only  if  the  difference 
between  the  wave  numbers  of  the  surface  waves,  5k,  is  equal  to 
the  wave-number  of  the  internal  wave.  Furthermore,  from  the 
second  equation,  If  5k  * 1 6k  | * It,  , 

. 8k  • U n, 

6 a ~ ~s  ^ 

6k  5k  * k3  . ... 


or,  in  terms  of  the  group  velocity  of  the  surface  waves, 

c (k)  + U cos  cp  = C • [5*21] 

g s 

where  cp  is  the  angle  between  the  surface  current  and  the  internal 
waves.  From  the  theory  of  resonant  interactions,  if  initially 
the  energy  resides  in  one  of  the  surface  wave  trains  (say  1)  and 
the  internal  v/ave  (3)*  the  amplitude  of  the  other  surface  wave 
train  (2)  at  a neighboring  wave  number  will  grow  — the  initially 
uniform  surface  wave  train  will  develop  modulations  or  groups 
that  grow  in  time.  This  is  the  phenomenon  of  resonance  that  has 
been  investigated  in  this  context  by  Lewis  et  al.  (1974)  for  the 
case  U **  0.  j 

It  is  important  to  realize  that  resonance  and  blockage  are 
two  different  ways  of  describing  the  same  physical  effect,  but 
each  description  has  its  own  limitations.  The  idea  of  resonance 
assumes  small  amplitude  waves  and  small  disturbances  — it  cannot 
cope  with  surface  velocity  variations  of  the  same  order  as  C, 
the  speed  of  the  internal  wave.  On  the  other  hand,  there  is  no 
restriction  on  the  relative  wave-  lengths  involved.  The  analysis 
of  Section  5.1  demands  that  the  surface  wave  length  be  short 
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compared  with  the  internal  wave  length  but  imposes  no  restric- 
tion on  the  magnitude  of  the  variation  in  U . It  is  not  intrin- 

s 

sically  a steady  state  approach.  When  (a)  the  surface  wave 
length  is  Indeed  much  less  than  the  internal  wave  length  (true 
under  oceanic  conditions)  and  (b)  the  variation  in  surface  cur- 
rent is  very  small  (sometimes  true,  sometimes  not  under  natural 
conditions),  the  two  approaches  overlap:  the  condition  [5.21] 

for  resonance  is  identical  to  [5-7]  for  blockage  (note  that  in 
[5.7]  only  the  case  = 0 is  considered). 

The  important  work  by  Lewis  et  al.  (1974)  has  (a)  provided 

a detailed  analysis  for  the  resonance  when  U = 0,  calculating 

s 

the  interaction  coefficients  and  (b)  verified  by  laboratory 
measurements  the  accuracy  of  the  theory,  showing  that  under 
resonant  conditions  the  modulations  in  energy  density  initially 
develop  such  that 


6E 

E 


[5.22] 


where  u is  the  maximum  orbital  velocity  of  the  internal  wave 
o 

and  k^  its  wave  number  and  finally  (c)  as  in  other  resonant  in- 
teraction experiments  (e.g.  McGoldrick,  Phillips,  Huang  and 
Hodgson,  1966)  the  width  of  the  resonance  band  decreases  as  the 
interaction  time  increases,  provided  the  internal  wave  slope  is 

small.  When  (U  - U )/C  is  not  negligibly  small,  the  "band 
s s 

width"  of  the  response  is  not  small  (see  Figure  12). 

5 . 4 Resonance  and  Blockage  in  the  Field 

The  work  of  Gargett  and  Hughes,  Phillips,  Lewis,  Lake  and 
Ko  has  established  a good  theoretical  foundation  for  these 
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phenomena  (If  not  a complete  one  in  all  respects)  and  the  ex- 
periments of  Lewis  et  al.  have  shown  that  the  resonance  phenom- 
enon can  be  reproduced  In  the  laboratory,  Gargett  and  Hughes 
described  some  suggestive  field  phenomena,,  but  qualitatively 
only,  Zonation  of  small  scale  waves  in  the  field  has  also  been 
established  under  some  conditions  (Perry  and  Sehimke,  1965): 
can  this  be  attributed  directly  to  resonance  or  blockage? 

Consider  first  of  all  the  wave  lengths  involved.  The  phase 
speed  of  internal  waves  in  the  ocean  Is  of  the  order  1 m/sec  at 
most;  the  estimate  of  Section  3*1  gives  0.7  m/sec.  The  wave 
length  In  meters  and  frequencies  In  cycles  per  second  (in 
brackets)  of  surface  waves  resonant  with  internal  waves  are 
given  in  the  table  below  for  various  values  of  the  surface 


current: 


TABLE  1 


C * 0.75 

0. 

,50 

0.25 

m/sec . 

u COS  4*  OB'  0.75 
s 

- 

-0.15 

(3*33) 

-0.60 

(1-7) 

0.50 

0.15  (3.33) 

- 

0.15 

(3.33) 

0.25 

0.60  (1.7) 

0.15 

(3.33) 

- 

0 

1*35  (1.1) 

0.60 

(1.7) 

0.15 

(3.33) 

-0.25 

2.40  (.83) 

1.35 

(1.1) 

0.60 

(1.7) 

0 

« 

0 

1 

3.75  (*67) 

2.40 

(-83) 

1.35 

(1.1) 

-0.75 

5*4  (.55) 

3.75 

.(•67) 

2.40 

(.83) 

(m/s ec) 

Wavelengths,  in  meter 

3,  and  frequi 

encies, 

bracketed.  In 

hz  of 

surface  waves  that  are  resonant  with  internal  waves  with  phase 
speed  C (m/sec.)  in  the  presence  of  surface  drift  IT  (m/ sec.)  at 
an  angle  < to  the  direction  of  motion  of  the  internal  waves.  The 
positive  direction  Is  taken  in  the  direction  of  motion  of  the 
internal  waves;  negative  wavelengths  indicate  surface  waves 
moving  in  the  opposite  direction. 
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Recall  also  that  for  short  surface  waves  their  amplitude  in  an 
actively  wind  generated  wave  field  is  determined  by  their  con- 
dition for  local  stability  — increased  energy  input  will  result 
not  in  Increased  amplitude  but  in  increased  wave  breaking,  (For 
wave  components  near  the  spectral  peak,  non-linear  conservative 
energy  exchanges  are  also  important  (Hasselmann  et  al.,  1973). 

An  energy  loss  will,  however,  result  in  a decreased  amplitude 
provided  the  local  rate  of  loss  of  energy  resulting  from  the 
interaction  exceeds  the  rate  of  energy  input  from  the  wind. 

If  these  short  wave  components  are  already  saturated  by  the 
wind  in  the  absence  of  an  internal  wave  disturbance,  under  what 
conditions  will  such  a disturbance  produce  modulations?  Clearly, 
when  the  wind  is  unable  to  regenerate  locally  the  wave  energy  in 
the  regions  where  it  is  decreasing  as  a result  of  the  resonant 
or  blockage  Interaction.  From  [5.22]  the  maximum  rate  of  de- 
crease of  wave  energy  is  given  by 


1 dE 
E dt 


-5Vi 


[5.23] 


where  k.  is  the  internal  wave  number  and  u is  the  magnitude  of 
i o 

the  surface  velocity  variations,  described  in  Section  3-  Under 
the  most  favorable  circumstances  for  resonance,  with  a mean  sur- 
face drift  nearly  equal  to  the  internal  wave  speed*,  uq  - C,  the 
speed  of  the  internal  wave,  so  that 


1 dE 
E dt 


< 


♦Resonance  theory  is  then  unreliable  for  such  large  surface 
currents. 


! 

? 
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where  n^  Is  the  radian  frequency  of  the  internal  wave.  A very 
high  frequency  internal  wave  on  the  thermocline  has  a period  of 
5 minutes  or  so  (this  corresponds  to  the  conditions  calculated  in 
Section  3.1)  so  that  n^  < 2 x 10_s.  Consequently,  for  the  most 
favorable  conditions  for  resonance  in  the  ocean, 

| ||  < -10'1  sec“l  [5.24] 

The  information  available  on  the  rate  of  energy  input  from 
the  wind  has  been  reviewed;  recent  relevant  measurements • have 
been  made  by  Wu  (1973).  Contours  of  the  growth  rate  of  surface 
waves  In  wind  are  shown  in  Figure  13.  For  surface  wave  lengths 
less  than  about  15  centimeters  that  influence  radar  backscatter- 
ing  measurements,  the  growth  rate  for  the  wind  exceeds  10"1  sec-1 
except  for  relatively  light  winds  (u#  < 25  cm/sec,  U < 5 m/sec 
— 10  knots);  if  the  wind  is  comparable  to,  or  greater  than  this 
limit,  modulations  in  the  short  wave  lengths  will  not  develop 
directly  from  resonant  interactions. 

Longer  wave  lengths,  whose  growth  rates  under  the  influence 
of  the  wind  are  slower,  may  develop  modulations,,  either  resonant 
or  non-resonant,  but  these  will  not  be  observable  directly  using 
surface  scanning  radar  with  wave  lengths  in  the  C,  X or  P bands. 
Nevertheless,  it  will  be  seen  in  Section  6 that  such  modulations 
may  In  turn  Influence  much  shorter  components  which  are  directly 
observable  in  this  way. 

Wave-wave  interactions  will  also  serve  to  make  more  uniform 
the  distribution  of  surface  wave  energy  in  the  face  of  Internal 
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wave  modifications.  These  have  been  considered  by  Zachariasen 
(1973)  under  the  assumption  that  energy  input  from  the  wind  and 
energy  dissipation  by  breaking  and  so  forth  are  linear  functions 
of  the  wave  energy  density.  It  is  difficult  to  draw  from  this 
analysis  any  simple  conclusions  except  that  (1)  if  the  inter- 
actions are  weak,  under  non-resonant  conditions  the  perturbations 
in  energy  density  are  in  phase  with  the  surface  current  as  in  a 
single  wave  train  (see  equation  5-12);  (2)  under  resonant  con- 
ditions there  will  be  perturbations  in  energy  density  in  phase 
with  the  rate  of  surface  strain  rather  than  the  surface  current, 
as  the  simple  (no  wind,  no  interaction)  theory  indicates.  The 
occurrence  of  these  modulations  in  the  theory  is,  however,  a 
consequence  of  the  "softness"  of  the  dissipation  function  assumed 
(linear  in  energy  density  E);  it  is  more  likely  that  the  rate  of 
dissipation  increases  more  rapidly  with  E than  linearly  and  this 
would  tend  to  suppress  the  modulations.  The  simple  equilibrium 
range  idea  is  essentially  equivalent  to  zero  dissipation  when 

E < E with  breaking  and  very  rapid  dissipation  when  E > E 

max  max 

•j 

5-5  Summary 

1.  Resonance  and  blockage  are  clearly  established  phenom- 
ena both  theoretically  and  experimentally,  in  the  laboratory. 

2.  In  the  field,  significant  resonance  effects  are  con- 
fined to  relatively  long  wave  lengths  (1-5  m or  so)  and  may 
occur,  the  optimum  conditions  being  (a)  relatively  large  internal 
wave  speeds  in  the  direction  of  surface  wave  travel  (b)  a sur- 
face current  opposed  to  the  direction  of  surface  wave  travel  and 
(c)  relatively  light  winds.  Such  modulations  will  not  be  di- 
rectly observable  using  conventional  surface-scanning  radar. 
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3.  Except  in  very  light  winds,  resonant  modulations  in 
short  waves  {wave  lengths  of  15  centimeters  or  less)  will  be 
masked  by  energy  input  from  the  wind,  and  any  modulations  that 
are  observed  by  radar  must  be  produced  by  some  other,  indirect 
effect,  in  which  resonance  may  play  a part,  but  not  the  sole 
part. 

4.  Amplitude  modulations  are  produced  by  non-resonant 
interactions  between  the  internal  wave  and  surface  waves  moving 
much  more  rapidly  than  the  Internal  wave.  Detailed  results  are 
given  in  Section  5-2.  Again,  these  modulations  are  not  directly 
observable  by  surface  scanning  radar  but,  as  will  be  seen  in 
Section  6,  can  produce  variations  in  small  scale  surface  wave 
structures  that  are. 
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6 . WAVE  BREAKING 

A breaking  wave  is  defined  simply  as  one  in  which  some  fluid 
elements  near  the  crest  are  moving  forward  faster  than  the  wave 
profile.  Wave  breaking  occurs  on  waves  over  a considerable  range 
of  scales;  at  the  larger  scales,  breaking  is  more  vigorous,  air 
is  entrained  and  the  breaking  is  clearly  visible  as  the  formation 
of  white  caps.  Small  scale  breaking  (wave  lengths  5-20  cm)  is 
however  widespread  in  an  active  wind-generated  wave  field  (Banner 
and  Phillips,.  1974);  a^r  n°4  usually  entrained,  the  flov;  being 
characterized  by  the  occurrence  of  surface • stagnation  points  and 
the  profile  by  the  appearance  of  irregular  'steps'  at  the  edge 
of  the  small  breaking  zone. 

6.1  Incipient  Freaking  • ’ . 


The  maximum  amplitude  that  a wave  train  (or  a local  group 
of  waves)  can  attain  without  breaking  depends  strongly  on  the 
wind  drift  in  the  upper  few  millimeters  below  the  air-water  inter- 
face. If  the  local  surface  drift  is  q (the  value  at  the  mean 
water  level  £ = 0;  very  closely  the  local  mean  value)  and  the 

phase  SDeed  of  the  wave  is  c,  then  the  maximum  amplitude  £ that 

' m 

the  wave  can  have  without  breaking  is  given  by  • ~ 


[6.1] 


(Banner  and  Phillips,  197*0,  as  illustrated  In  Figure  14  . In 
the  absence  cf  surface  drift,  f = c2/l;g  the  Stokes  vaxue,  but 
under  natural  conditions  the  maximum  wave  amplitude  is  substan- 
tially less,  'lie  effect  is  particularly  significant  at  those 
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short  wavelengths  that  are  responsible  for  radar  backscattering. 
In  a wind  of  10  m/sec , q^  **  30  cm/sec.  A surface  v;ave  length  of 
10  centimeters  has  e * 40  cm/sec  and  s maximum  amplitude  of  only 
0.03  times  the  Stokes  limit.  If  the  locally  averaged  value  of 
the  surface  drift  varies  by  an  amount  5a,  the  maximum  wave  ampli- 
tude varies  as 


C 


m 


[6.2] 


so  that  for  wave  speeds  only  slightly  greater  than  q , even  small 

increases  in  q will  produce  substantial  reductions  in  r . 

m 

If  the  local  wave  amplitude  is  forced  to  exceed  the  limit' 
[6.1],  water  at  the  crest  tumbles  forward,  the  wave  breaks,  losing 
energy  until  the  local  energy  density  is  sufficiently  reduced 
that  breaking  ceases.  It  is  suspected,  but  not  known  with  cer- 
tainty, that  there  is  a hysteresis  involved  in  breaking;  waves 
that  begin  to  break  may  continue  to- do  so  until  their  amplitude 
is  significantly  less  than  [6.1].  The  onset  of  breaking  in  a 
particular  wave  crest  is  generally  quite  abrupt;  there  is  a non- 
linear feedback  process  resulting  from  the  compression  of  the 
vortical  layer  that  rapidly  augments  the  forward  surface  surge 
as  the  breaking  condition  is  approached. 


Once  waves  on  a small  scale  begin  to  break,  as  in  Figure  16, 
the  wave  profile  is  far  from  sinusoidal.  The  breaking  region 
ahead  of  the  crest  has  an  irregular  area  of  high  slope  where  the 
water  is  falling  forward;  in  terms  of  the  Fourier  components, 
many  higher  harmonics  are  generated.  The  high  frequency  unsteady 
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motion  in  the  small  tumbling  region  also  generates  many  capillary 
components  with  wave  lengths  of  a few  millimeters  that  may  not 
be  excited  directly  by  the  wind. 

The  presence  of  a surface  slick  that  is  saturated  or  close 
to  saturated  has  a very  strong  influence  on  wave  breaking.  Com- 
pression of  the  slick  is  resisted  by  the  gradient  of  the  film 
pressure  (see  Section  4),  so  that  the  forward  surge  near  the 
short  wave  rrast  that  leads  to  breaking  is  resisted  and  breaking 
is  delayed.  Energy  dissipation  is,  of  course,  augmented  even 
without  breaking,  though  it  seems  to  be  the  inhibition  of  short 
wave  breaking  that  gives  the  dramatic  difference  in  surface 
appearance  between  a concentrated  slick  and  a surrounding  cleaner 
surface.  No  longer  are  small  scale  capillary  wavelets  generated 
by  the  breaking  and  the  attenuation  rate  is  generally  sufficient 
to  prevent  direct  generation  from  the  wind. 

6 .2  The  Suppression  of  Chert  Waves  by  fwell  in  the  Presence  of 

Wind  Drift 

The  surface  wind  drift  is  not  uniform  in  the  presence  of  a 
wave  field  because  of  the  surface  straining  associated  with  the 
waves.  When  waves  are  propagating  generally  in  the  same  direc- 
tion as  the  drift,  the  drift  is  augmented  in  the  vicinity  of  the 
wave  crest.  The  wave  breaks  when  the  augmentation  is  such  that 
the  surface  speed  at  the  wave  crest  exceeds  the  phase  speed  of 
the  wave,  but  the  amplification  is  present  even  in  non-breaking 
waves.  The  maximum  value  of  the  surface  drift  at  a long  wave 
crest  is  given  by  Phillips  and  Banner  (1974)  as 


HYDROHAUTlCS,  Incorporated 

-56- 

« * <Ct  * u>  - {(ct  • “Is  - ’o  - %>}*  [6-3] 

where  C.  is  the  long  wave  speed,  u its  orbital  speed  at  the  sur- 
face  (=  C ak  where  a * long  wave  amplitude,  k * wave-number)  and 

Kt 

is  the  value  of  the  surface  drift  at  points  where  £ * 0 on  the 
mean  water  level,  or  approximately  the  average  value  of  the  sur- 
face drift  over  the  whole  surface.  For  wave  periods  greater  than 
about  2 seconds  (X  > 10  m)  variations  in  wind  stress  may  influence 
the  distribution  (see  Section  2 of  this  report) 5 the  effect  should 
be  calculated  numerically  but  it  is  expected  that  the  Increased 
wind  stress  at  the  wave  crests  may  even  augment  further  the  in- 
creased drift  there  produced  by  the  dynamics  of  the  strained 
surface  layer.  The  neglect  of  this  further  effect  should  repre- 
sent an  under-estimate  of  the  effects  on  short  waves. 

When  short  waves  are  being  overtaken  by  a longer  one,  near 
the  long  wave  crest  the  short  waves  experience  an  augmented  wind 
drift  qQ,  so  that  the  maximum  height  that  they  can  have  at  in- 
cipient breaking  is  reduced  in  accordance  with  equation  [6.1). 

In  fact,  in  an  active  wind  with  continuing  energy  supplied  to 
the  short  waves,  they  can  be  expected  to  actually  break  on  a 
small  scale  as  they  approach  the  long  wave  crest  and  the  local 
mean  drift  increases  so  that  [6.1]  becomes  more  stringent. 

Once  the  long  wave  crest  has  passed  they  may  be  regenerated  by 
the  wind  but  only  to  break  at  the  next  long  wave  crest  when  the 
Increased  drift  again  limits  the  amplitude.  Accordingly,  the 
presence  of  longer  waves  and  wind  drift  reduce."  the  amplitude 


HYDROuAUTICS,  Incorporated 

-57- 

of  short  waves,  the  amount  of  the  reduction  being  dependent  on 
the  long  wave  speed  and  slope.  By  assuming  that  the  average 
energy  density  in  a short  wave  field,  constantly  regenerated  by 
the  wind,  is  proportional  to  (in  fact,  roughly  equal  to)  the 
energy  density  when  the  short  waves  are  at  the  point  of  incipient 
breaking  at  the  long  wave  crest,  Phillips  and  Banner  (197^+)  were 
able  to  account  quantitatively  for  reductions  in  short  wave  energy 
density  produced  by  long  waves  in  a variety  of  experimental  sit- 
uations, Early  measurements  by  Mitsuyasu  (1966)  were  also  very 
consistent  with  the  theoretical  predictions  (see  Figure  17) 
though  experimental  verification  is  still  rather  sparse. 

For  short  waves  superimposed  on  long  ones,  the  maximum 
amplitude  at  incipient  breaking  Is  given  by 

<Bax  * <cc  * q)S  [6'4] 

where  q represents  the  locally  averaged  surface  drift  [6.3]  in- 
duced by  the  long  wave  at  its  crest.  This  can  be  expressed 
alternatively  as 

I 

A , (1  . ak)  - ju  - ak)2  - “ (a  - T?]}  [6-5] 

where  qQ  is  the  overall  average  surface  drift  and  ak  the  long 
wave  slope,  c is  the  short  wave  speed  at  the  long  wave  crest 
(also  modified  by  the  straining  of  the  long  wave  from  its  un- 
disturbed value  c ): 

o' 

c (1  - ak) 


c « 
c 


[6.6] 
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and  gf  is  the  apparent  gravitational  acceleration  experienced  by 
the  short  wave 

• - — - . g»  ss  g(i  . ak) . [6.7] 

6.3  Short  Wave  Modulations  Associated  with  Internal  Waves 

J 

Assume  that  internal  waves  produce  a distribution  of  surface 

current  whose  variations  may  or  may  not  be  augmented  by  a ;mean 

* 1 

surface  current  as  in  Section  3 of  this  report.  Assume  'that  the 
wind  is  blowing  in  some  direction  generating  waves  over  a range 
of  scales;  let  us  (rather  simplistically ) represent  these  waves 
initially  as  a train  of  long  waves  moving  much 'faster  than  the 
internal  wave  superimposed  on  which  are  short  wavelets  (say  in 
the  2-10  cm  range).  The  characteristic  growth  time  for  the  long 
waves  is  taken  to  be  large  (much  larger  than  the  time  to  propa- 
gate from  one  internal  wave  crest  to  the  next)  so  that  the  direct 
effect  of  wind  re-generation  can  be  neglected  during  these  time 
intervals.  On  the  other  hand,  the  growth  time  for  the  short 
waves  is  small;  they  are  generally  saturated  and  if  thei^  ampli- 

, l 1 

tude  is  severely  restricted  at  a long  wave  crest,  they  will 
rapidly  be  regenerated  once  the  crest  is  passed.  The  quantita- 
tive form  of  these  assumptions  is  considered  later. 

Since  the  response  time  of  the  wind-drift  layer  (in  the 
upper  few  milimeters)  is  short  (see  Section  2)  the  mean  wind 
drift,  averaged  over  an  area  small  compared  with  an  internal 
wave  length  but  large  compared  with  a long  surface  wave  length 
is,  in  essence,  constant.  The  direct  straining  of  the  surface 
wind-drift  layer  by  the  internal  wave  is  too  slow  and  of  too 
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small  a magnitude  (by  several  orders)  to  modify  the  structure  of 
the  wind  drift  layer  directly.  The  short  time  scales  appropriate 
to  this  layer  Indicate  that,  so  far  as  the  internal  waves  are 
concerned,  the  surface  drift  layer  is  everywhere  in  local  equi- 
librium with  the  wind  and  shorter  waves. 

The  long  waves  propagate  relatively  rapidly  through  the 
variable  current  pattern  produced  by  the  internal  waves.  They 
are  not  resonant  and  they  are  not  blocked,  but  (relatively  small) 
modulations  in  long  wave  amplitude  are  produced  as  described  in 
Section  5-  The  long  waves  are  smallest  when  the  surface  drift 
component  in  the  direction  of  internal  wave  travel  is  largest 
(see  Eq.  5.18}  if  0 < tt/2.  The  small  modulations  In  long  wave 
slope  produce  modulations  in  the  maximum  values  of  surface  drift 
at  the  long  wave  crests,  where  it  is  augmented  by  the  convergence 
in  the  long  waves.  An  important  point  is  that  the  surface  drift 
generated  at  the  long  wave  crest  is  dependent  on,  and  a substan- 
tial fraction  of  the  long  wave  speed,  so  that  the  modulations 
in  surface  speed  at  the  different  long  wave  crests  are  a (fairly 
small)  fraction  of  the  long  wave  speed.  However,  they  are  a 
substantial  fraction  of  the  short  wave  speed . Short  waves  at 
different  phases  of  the  internal  wave  experience  surface  drifts 
at  the  long  wave  crests  that  vary  by  a substantial  fraction  of 
their  own  speed;  in  the  previous  section  it  was  shown  that  the 
maximum  amplitude  the  short  waves  can  have  against  breaking  and 
in  the  presence  of  a regenerating  wind,  depends  sensitively  on 
the  value  of  the  drift  q as  a fraction  of  their  phase  speed  c. 
Consequently,  there  will  be  modulations  in  the  short  wave  energy 
density,  with  maxima  where  the  long  waves  are  smallest,  that  is. 
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when  the  surface  current  in  the  direction  of  wave  propagation 
is  greatest.  Minima  will  occur  when  the  current  is  least  or  in 
the  reverse  direction. 

These  statements  will  now  he  put  into  more  quantitative  form 
The  variation  In  (root  mean  square)  amplitude  of  the  long  waves  • 
(no  resonance)  is  given  by  [5.17]  as 


5a  s j 1 + cos  2g  - y(2  + cos  e) 

a “ ” ^ cJ  (cos  6 - 2y ) (1  - 2y) 

and  the  variation  in  wave-number  by  [5.14]  as 


2 cos  $ 
1 - 2y 


[6.9] 


where  y « -(Ug  - C)/C^,  Uc  being  the  component  of  the  mean  sur- 
face current  (in  the  top  few  meters)  in  the  direction  of  the  long 
waves.  The  slope  variation  is  accordingly 


5 (ak) 
ak 


4 cos8  9 - y(2  4-  5 cos  o) 
(cos  e - 2y ) (1  - 2y) 


[6.10] 


Modulations  in  long  wave  slope  when  SU^/C  * 0.2  are  shown  In 
Figure  13  as  a function  of  C /C.  The  resonance  point  is  when 

T«« 

C /C  * 2;  we  are  concerned  with  considerably  longer  waves  such 
that  C./C  > 2.  The  slope  modulations,  opposed  in  phase  to  the 
surface  current,  decrease  as  C,/C  increases;  a characteristic 

V 

magnitude  is  0.2  or  20£. 

Now  from  [6.5],  taking  differential  Increments,  the  modula- 
tions in  surface  wind  drift  at  the  crests  of  the  long  waves  are 
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given  as  a fraction  of  the  long  wave  speed  by 

ft  q q \1*"2 

P - e(ak)  / (1  - ak)2  - ^ 2 - ^ (1  - ak)  - 1 

■t  I K I 


Note  that  qQ,  the  average  surface  wind  drift  is  essentially  con- 
stant with  respect  to  phase  of  the  internal  wave  on  account  of 
the  short  response  time  of  the  wind  drift  layer  — see  Section  2. 
Note  also  that  in  the  field,  q /C  « 1,  the  average  surface 

O 

wind  drift  is  snail  compared  with  the  long  wave  speed.  The  ex- 
pression [6.11]  is  illustrated  in  Figure  19  for  various  values 
Of  d0/Ct. 

Consider  now  the  short  wavelets,  constantly  regenerated  by 
the  wind,  riding  on  these  long  ones.  The  maximum  amplitude  they 
can  have  against  breaking  is  given  by  [6.4],  i.e. 


«n,ax  * (o  - q)s 


where  g'  * g(l  - ak)  is  the  apparent  gravity  acting  on  the  short 
waves  as  they  ride  over  the  long  wave  crest,  c * cq  (1  - ak)  is 
the  local  wave  speed  of  the  short  waves  and  q is  the  local  sur- 
face wind  drift  at  the  long  wave  crest.  Again  taking  differ- 
ential increments,  the  variations  are  given  by 
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*max 


cc-T-a 


s’ 

+ c-a  1-ak 
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1-ak 
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1 
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o 

1 
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2 f (ak)  c 


1-ak  (c-q)(l-ak) 


\ 6a 
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[6.12] 
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In  this  expression,  £(ak)  and  5q  are  both  directly  proportional 

to  5U  according  to  [6.11]  and  [6.10];  in  the  calculations  below, 
s 

6U  is  taken  as  0.2C  (a  rather  low  value;  see  Figure  5)  so  that 
s 

for  other  values  of  the  ratio,  the  modulations  can  be  found  by 
simple  multiplication. 

Since  the  modulations  In  short  wave  energy  density  5E/E  * 
8(6tm)/c«x  they  are  «Ivcn  iy 


&E 

E 


-2ak  {Fig.  16] 


c + q 

o 


(1-ak)  (co-q) 


£ c. 


vq 


{Fig.  19} 


[6.13] 


where  {Fig.  18],  {Fig.  19}  represent  the  ordinates  for  appro- 
priate conditions  from  these  curves.  Note  that  when  q - cq,  the 
amplitude  of  the  proportional  modulations  becomes  indefinitely 
large,  the  short  waves  being  erased  from  at  least  part  of  the 
longer  wave  cycle. 

The  results  of  some  detailed  calculations  based  on  these 
expressions  are  shown  in  Figures  20  through  27.  In  interpreting 
these  figures,  several  points  should  be  kept  in  mind. 


1.  In  each  case,  the  ordinate  represents  the  amplitude  of 
the  modulations  In  .short  wave  energy  density  over  a period  of  the 
internal  wave. 


2.  The  short  waves  considered  here  are  freely  travelling 
waves  under  the  influence  of  wind,  not  the  Fourier  harmonics  of 
longer  waves  whose  basic  wave  length  is  a multiple  of  the  wave 
length  given.  For  freely  travelling  short  waves  of  a given  wave 
length,  there  Is  a cut-off  produced  by  the  wind  drift,  while 
longer  waves,  with  harmonics  of  the  same  wave-number  can  still 
propagate. 
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3.  Calculated  variations  of  AE/E  larger  than  some  limit  of 

order  1 can  no  longer  be  regarded  as  perturbations  and  the  num- 
erical calculations  are  inaccurate.  Indicated  modulations  great- 
er than  this  simply  indicate  a very  large  response. j Values  of 
AE/E  less  than  10 % may  possibly  be  regarded  as  insignificant  and 
are  not  plotted.  | 

4.  All  calculations  $re  for  surface  current  variations 

1 ! 

whose  amplitude  is  0.2  tim^i  the  internal  wave  speed.  For  other 

i 

values,  the  curves  should  be  increased  or  decreased ( in  linear 
proportion. 

5.  Under  natural  conditions,  the  wind  wave  spectrum  is 
generally  narrow  and  sharply  peaked.  Ine  complex  surface  wave 
field  is  here  simplified  to  a long  wave,  together  with  super- 
imposed local  short  waves  limited  by  breaking. 

» 

> Figure  20,  21  and  22  give  AE/'E  for  2.5,  10  and  20  centimeter 
waves  as  a function  of  wind  speed  for  the  case  of  no  mean  surface 
current,  long  wave  slope  0.(2,  for  various  values  of  the  wave 
length  of  the  longest  waves!,  * The  bracketed  pairs  of  curves  show 
the  influence  of  variations  in  internal  wave  speed  from  30  tc  J:0 
cm/sec.  Note  the  cut-off  in  these  curves  when  the  surface  win  1 
drift  is  sufficient  to  suppress  freely  travelling  waves  (qQ  = c). 

Figures  23,  24,  25  give  AF/E  again  for  2.5,  10  and  20  centi- 
meter waves  as  a function  of  the  slope  of  the  longest  waves, 
again  with  no  mean  surface  current  and  with  the  wind ' speed  of 
7-5  m/sec  (15  knots).  ?,Tote  the  much  greater  response  of  shorter 


wave  s . 
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Figure  2 6 gives  $E/E  for  the  three  short  wave  lengths  as  a 
fuction  of  mean  surface  current  referred  to  the  internal  wave 
speed  (though  with  a constant  value  of  0,2  for  the  variations  in 
surface  current  to  internal  wave  speed,  that  are  produced  by  the 
internal  waves).  The  longer  wave  has  1 m wavelength,  internal 
wave  speed  20  cm/sec  and  again  a 7.5  "m/sec  (15  knot)  wind. 

Figure  27  gives  the  directional  dependence  of  the  response 
as  a function  of  angle  0 between  the  long  surges  waves  and  the 
internal  wave.  There  is  some  directional  weakening,  but  the 
response  is  surprisingly  broad  band. 

The  general  trends  of  the  calculations  are  Indicated  below. 
6.4  Suxma  ry 

1.  The  response  of  2.5  cm  waves  is  larger  than  that  of  10 
or  20  cm  waves  by  a factor  of  order  3 (at  lov;  wind  speeds, 

U — 2 m/sec,  or  ^ knots)  to  10  or  much  more  at  higher  wind  speeds 
(greater  than  7-5  m/sec,  15  knots). 

2.  The  phase  of  the  short  ..  ve  minimum  corresponds  to  the 
internal  wave  phase  where  the  current  has  its  maximum  adverse 
value  with  respect  to  the  longer  waves  in  the  surface  wave  field 
(above  the  crests  if  the  two  are  concurrent).  The  maximum  occurs 
between  this  point  and  the  following  internal  wave  trough. 

3.  Orders  of  magnitude  of  variations  with  the  various 
parameters  are  as  follows; 
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7 . EPILOGUE 

This  work  has  sought  to  evaluate  the  predominant  influences 
of  internal  waves  on  short  surface  waves  under  natural  condi- 
tions where  a variety  of  physical  processes  occur  In  combination 
The  strongest  and  most  widespread  effect  appears  to  arise  from  a 
relatively  Small  modulation  to  long  waves  that  run  through  the 
internal  wave-induced  motion,  these  modulations  having  a some- 
times very  strong  influence  on  the  maximum  height  that  short 
waves  can  maintain  against  energy  loss  by  breaking  while  under 
the  continued  Influence  of  the  wind.  The  results  are  summarized 
in  Section  6.4.  Modulations  in  short  wave  energy  density  and 
in  density  of  short  wave  breaking  and  surface  stirring  are  im- 
portant. Since  the  present  report  is  already  long  enough,  these 
questions  will  be  considered  in  a sequel. 
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FIGURE  1 - THE  VELOCITY  GRADIENT  AT  THE  WATER  SURFACE  AS  A FUNCTION 
OF  THE  SQUARE  OF  THE  CENTERLINE  SPEED.  The  linearity  confirms 
that  the  shear  stress  Is  supported  predominantly  by  molecular  viscosity. 
(After  Wu,  1968) 


FIGURE  2 - SURFACE  DRIFT  CURRENT  AS  A FUNCTION  OF  CENTERLINE  WIND  SPEED 
Measurements  by  Wu  (1968)  and  Banner  and  Phillips  (1974)  in  two  different 
wind-wave  facilities. 
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FIGURE  5 - DISTRIBUTIONS  OF  SURFACE  CURRENT,  WITH  RESPECT  TO  PHASE 
OF  AN  INTERNAL  WAVE,  WITH  o/d  =0.1.  r = U /C,  the  ratio  of 
mean  current  to  internal  wove  speed.  Note  that  the  magnitude  of  the 
modulations  increases  with  r;  above  r = 0.5641  in  this  cose,  fluid 
above  the  wave  trough  moves  faster  than  the  wave  profile  and  a reg- 
ion of  closed  mean  circulation  develops.  r = 0 corresponds  to  the 
classic  irrotational  case. 
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FIGURE  7 - THE  AMPLIFICATION  IN  SURFACE  CURRENT  MODULATIONS 
RESULTING  FROM  THE  EXISTENCE  OF  A MEAN  CURRENT. 
Above  the  broken  line,  fluid  at  the  surface  above  the  internal 
wave  troughs  moves  forward  faster  than  the  v/ave  itself  and  a 
closed  eddy  forms. 
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FIGURE  8 - CONTO  IRS  OF  THE. AMPLIFICATION  FACTOR  OF  SURFACE  CURRENT 
MODULATIONS  (solid  curves, identified  by  numbers  in  circles)  AS  A 
FUNCTION  OF  a/d  AND  0S/C  AND  OF  THE  TOTAL  AMPLITUDE  OF 
CURRENT  MODULATION  DIVIDED  BY  INTERNAL  WAVE  SPEED  (broken 
curves, identified  by  numbers  in  boxes). 
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FIGURE  11  - THE  CONDITIONS  FOR  FORMATION  OF  A SATURATED 

SLICK  IN  GENERALLY  DISPERSED  SURFACE  MATERIAL.  Slicks 
will  be  formed  under  conditions  that  correspond  to  point* 
inside  the  bowl  if  Vw  + q = 0 , or  inside  tbe  V -shaped 
region  if  tbe  surface  current  vanishes,  o - internal  wove 
amplitude,  d = tbermocline  depth. 
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FIGURE  12  - BLOCKAGE  CONDITIONS  FOR  SURFACE  WAVES  OF  SPEED  c WITH 
INTERNAL  WAVES  CF  SPEED  C AND  A SURFACE  DRIFT,  WHOSE 
MEAN  VALUE  IS  U$  AND  WHOSE  MINIMUM  IN  THE  PRESENCE  OF 
OF  INTERNAL  WAVES  IS  Lls  . Note  that  when  the  variations  in  sur- 
face ^current  Us  — Us  - 0,  blockage  is  reduced  to  a single  point: 

C ~ Us  = (?)  cQ  . the  resonance  condition  (5.8). 
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FIGURE  15  - STREAK  LINES  IN  A SMALL  SCALE  BREAKING  WAVE 
For  Identification,  see  Figure  16. 
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FIGURE  16  - AN  IDENTIFICATION  OF  FIGU 
The  undisturbed  water  depth  is  1! 
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FIGURE  17  - THE  REDUCTION  IN  ENERGY  DENSITY  OF  SMALL-SCALE 
BREAKING  WAVES  PRODUCED  BY  SUPERIMPOSED  LONGER 
WAVES  AFTER  PHILLIPS  AND  BANNER  (1974).  The  measure- 
ments are  by  Mitsuyasu  (1966);  the  line  represents  the  theoret- 
ical calculation,  with  no  adjustable  constants. 
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FIGURE  20  - PROPORTIONAL  MODULATIONS  IN  ENERGY  DENSITY  OF  SH 
BREAKING  WAVES  ASA  FUNCTION  OF  WIND  SPEED.  Short  w, 
length  * 2.5  cm,  surface  current  = 0,  long  wave  slope  = 0,2.  Th< 
four  pairs  of  curves,  bracketed,  are  for  long  wavelengths  of  1 me 
3 meters,  10  meters,  and  30  meters,  the  left-hand  member  of  eac 
pair  being  for  an  internal  wave  speed  of  40  cm/sec,  the  other  for 
30  cm/sec. 
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FIGURE  21  - PROPORTIONAL  MODULATIONS  IN  ENERGY  DENSITY  OF  SHORT 
BREAKING  WAVES  ASA  FUNCTION  OF  WIND  STEED.  Short  wave- 
length = 10  cm  , surface  current  = 0,  long  wove  slope  = 0.2.  The 
four  pairs  of  curves,  bracketed,  are  for  long  wavelengths  of  1 meter, 
3 meters,  10  meters,  and  30  meters,  the  left-hand  member  of  each 
pair  being  for  an  internal  wave  speed  of  40  cm/sec,  the  other  for 
30  cm/sec . 
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FIGURE  22  - PROPORTIONAL  MODULATIONS  IN  ENERGY  DENSITY  OF  SHORT 
BREAKING  WAVES  ASA  FUNCTION  OF  WIND  SPEED,  Short  wave- 
length = 20  cm  , surface  current  - 0,  long  wave  slope  = 0,2.  The 
four  pairs  of  curves,  bracketed,  are  for  long  wavelengths  of  1 meter, 
3 meters,  10  meters,  and  30  meters,  the  left-hand  member  of  each 
pair  being  for  an  Internal  wave  speed  of  40  cm/sec,  the  other  for 
30  cm/sec. 
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THE  FRACTIONAL  MODULATIONS  IN  ENERGY  DENSITY  FOR 
SHORT  WAVELENGTHS  OF  2.5,  10,  and  20  cm  AS  A FUNCTION 
OF  MEAN  SURFACE  CURRENT,  BUT  WITH  A CONSTANT  VALUE 
OF  0.2  FOR  THE  VARIATIONS  IN  SURFACE  CURRENT.  Longer 
wavelengths  = 1 m,  Internal  wave  speed  = 30  cm/sec,  wind  speed= 
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FIGURE  27  - THE  DIRECTIONAL  RESPONSE  IN  SHORT  WAVE  ENERGY  MODULATIONS 
AS  A FUNCTION  OF  ANGLE  BETWEEN  LONG  SURFACE  WAVES  AND  THE 
INTERNAL  WAVE.  When  this  angle  is  as  large  as  60°,  the  magnitude  of  the 
modulations  is  reduced  only  by  about  one-half. 


